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ABSTRACT: 

Two-particle momentum correlations between pairs of identical particles produced in relativistic 
heavy-ion reactions can be analyzed to extract the space-time structure of the collision fireball. 
We review recent progress in the application of this method, based on newly developed theo- 
retical tools and new high-quality data from heavy-ion collision experiments. Implications for 
our understanding of the collision dynamics and for the search for the quark-gluon plasma are 
discussed. 
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1 INTRODUCTION AND OVERVIEW 

1.1 Intensity Interferometry 

The method of two-particle intensity interferometry was discovered in the early 
1950's by Hanbury Brown and Twiss (HBT) |jl] who applied it to the measure- 
ment of the angular diameter of stars and other astronomical objects. These first 
measurements used two-photon correlations. An independent development oc- 
curred in the field of particle physics in 1960 by Goldhaber, Goldhaber, Lee and 
Pais Q who extracted from two-pion correlations the spatial extent of the annihi- 
lation fireball in proton-antiproton reactions. The method exploits the fact that 
identical particles which sit nearby in phase-space experience quantum statistical 
effects resulting from the (anti)symmetrization of the multiparticle wave function. 
For bosons, therefore, the two-particle coincidence rate shows an enhancement at 
small momentum difference between the particles. The momentum range of this 
enhancement can be related to the size of the particle source in coordinate space. 
HBT interferometry differs from ordinary amplitude interferometry that it 
compares intensities rather than amplitudes at different points. It shows the 
effects of Bose or Fermi statistics even if the phase of the (light or matter) wave 
is disturbed by uncontrollable random fiuctuations (as is, for example, the case 
for starlight propagating through the earth's atmosphere) or if the counting rate 
is very low. To illustrate how the method works and how its applications in astro- 
nomy and in particle physics differ let us consider the following simple model |P,Q|: 
two random point sources a and 6 on a distant emitter, separated by the distance 
R, emit identical particles with identical energies Ep = (m? -l-p^)"*^" which, after 
travelling a distance L, are measured by two detectors 1 and 2, separated by the 
distance d (see Figure Ilia). L should to be much larger than R or d. The total 
amplitude measured at detector 1 is then 

Ai = y (a e*(P'-i-+'^-) + /3 e^(P^»+'^'')) , (1) 

where a, (3 are the amplitudes emitted by points a and h, (pa, <Pb are their random 
phases, ria,ri6 are their distances to detector 1. The total intensity in 1 is 

h = ^ (\a\^ + |/?|2 + 2 Re a* 13 e*W'^"-^i")+(<^'>-'^")l) , (2) 

with a similar result for l2- The last term in (pi) contains information on the 
distance R between the sources a and 6, but it vanishes after averaging the signal 
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(a) (b) (c) 






Figure 1: Measurement of the separation R of two sources a and b by correlating 
the intensities in detectors 1 and 2 at varying distances d. a: The general scheme. 
b: The specific situation in astronomy, c: The specific situation in particle 
physics. 



over some time, i.e. over the random phases (j)a,b 

{Il) = {l2)- ^ "-'' 
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(3) 



The product of the averaged intensities (/i)(/2) is thus independent of both R 
and d. 

The same is not true for the time-averaged coincidence rate which is obtained 
by multiplying the two intensities before averaging: 



(/i/2) = (/i)(/2) + -^|a| 



cos 



[pina 



r2a - rib + r2b 



The two-particle intensity correlation function is thus given by 

(/1/2) , , 2|a|2|/3|2 



C{R,d) 



1 + 



■ cos 



[pina 



r2a - rib + r2b 



(/i)(/2) ^ ' (|aP + |/3|2)2 
For large L ^ R,d, the argument of the second, oscillating term becomes 

ria - r2a - rib + r2b — ' —j- (cos{d, R) - cos(d, L) cos{R, L)j . 



(4) 



(5) 



(6) 



Note the symmetry of this expression in d and R, the separations of the detectors 
and of the emitters; this symmetry is lost in the two practically relevant limits: 

1. In astronomical applications the emission points a and b are part of a star's 
surface or of even larger objects, while the detectors on earth are only a few 
meters or kilometers apart: R^ d. In this limit (see Figure ^3) the cosine- 
term in (^) reduces to cos{d-{pa—pb)), with Pa^b = P^a,b, the unit vectors 
^a,b giving the directions from the detectors to the two emission points a, b. 
Experimentally one varies the distance d between the detectors, and from 
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the resulting oscillation of the signal one extracts the angular separation 9 of 
the two emitters via \Pa—Pb\ ~ ^/-^- Absolute determination of the source 
separation R is possible only if their distance L is known from elsewhere. 

2. In nuclear and particle physics the emitters are very much smaller than 
the detector separation, i? <C d. Then (see Figure ||c) the cosine-term 
in (|5|) becomes cos(-R-(pi— P2))) in similar notation as above. Now the 
experimental control variable is the momentum difference Pi—p2, and R is 
extracted from the oscillation period. 

In real life one has, instead of two discrete ones, a continuum of sources described 
by a distribution p{R) of their relative distances. For the case R <ti d, after 
averaging (g) over this relative distance distribution, the measured correlation 
function is then given by the Fourier transform of p{R): 

C{pi -P2)-1^ Jd'R p{R) cos(i^ ■{P1-P2)). (7) 

As we will see, this expression is only applicable to static sources. The assumption 
of a static source is adequate for stars. The particles emitted in high energy 
hadron or nuclear reactions, however, come from sources which feature strong 
dynamical evolution both along the beam direction and transverse to it. As 
a result, two-particle correlation measurements in heavy- ion physics exhibit a 
much richer structure, but their interpretation is also considerably more involved. 
The present review covers the technical tools required to access this richness of 
information and their practical application in heavy- ion reactions. 

1.2 HBT Interferometry for Heavy-Ion Collisions 

High energy heavy-ion collisions produce hadronic matter under intense condi- 
tions of temperature and density. While the highest densities are reached in the 
early stages of the collision, most of the observed particles are hadrons, which 
are emitted rather late in the evolution of the collision. For this reason, the 
measured inomentum spectra and correlations contain direct information only 
about the size, shape, and dynamics of the source at "freeze-out", i.e. when the 
hadrons cease to interact. 

The dynamical information is of particular importance as it allows us to connect 
the observed final state with the hot and dense early stages. Therefore much of 
the effort in the last few years has gone into the extraction of this dynamics. 
It turns out that both the single-particle spectra and two-particle correlations 
are sensitive to certain combinations of thermal and collective motion in the 
source. A simultaneous analysis of particles with different masses allows for 
a separation of the two effects: while the thermal motion generates a common 
momentum distribution in the source, the collective motion creates a flow velocity 
which combines with the thermal momentum proportionally to the particle's 
mass. Further discrimination is achieved by combining the spectra with the 
two-particle correlations which reflect the collective source dynamics through a 
characteristic momentum dependence of their width parameters. 

The aim of HBT (or two-particle Bose-Einstein) interferometry is therefore to 
provide, in conjunction with a simultaneous analysis of single-particle spectra, a 
complete and quantitative characterization (both geometrically and dynamically) 
of the emitting source at "freeze-out" . Such a characterization can be used for 
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backward extrapolations into the interesting dense earlier stages of the collision, 
and it provides stringent constraints on dynamical models aiming to describe the 
kinetic evolution of these stages. 

1.3 New Developments During the Last Decade 

The last decade has brought great strides in measurement, theory, and inter- 
pretation of two-particle correlations. Dedicated experiments, optimized for mo- 
mentum resolution, allow measurement of correlations to very small momentum 
difference. Much increased statistics of particle pairs has opened the possibility 
of multidimensional correlation analysis. Correlation functions of identified kaons 
have become available; as many fewer kaons arise from resonance decays, they 
provide a more direct picture of the emitting source than the more prevalent 
pions. Correlations of proton pairs, which will not be discussed here, are also 
measured in heavy-ion collisions. The experimental programs at the CERN SPS 
and Brookhaven AGS have allowed comparison between small and large systems 
at different energies, using S and Pb beams at CERN and lower energy Si and 
Au beams at the AGS. The comparisons are greatly aided by the development of 
a commonly accepted analysis formalism among the experiments. Furthermore, 
the pair momentum dependence of the correlations is now being used to provide 
dynamical information about the space-time structure of the particle source. 

Intensive modeling with event generators, combined with methods to extract 
correlation functions from them, has been used to study experimental effects such 
as acceptance and to verify the usefulness of simple model parameterizations in 
theory. New parametrizations of the two-particle correlation function have been 
developed which are particularly well adapted for the sources created in high 
energy collisions. Recently, the necessity of separating geometrical, temporal 
and dynamical aspects of the correlations has been recognized, and methods to 
do so have been developed. Intensity interferometry has thus developed into a 
quantitative tool which at the present stage yields important source parameters 
at the level of 20% accuracy. 

2 THEORETICAL TOOLS FOR ANALYZING TWO-PARTI- 
CLE CORRELATIONS 

2.1 1- and 2-Particle Spectra and Their Relation to the Emitting Source 

The 2-particle correlation function is defined as the ratio of the Lorentz-invariant 
2-particle coincidence cross section and the product of the two single particle 
distributions: 

^, . ^ EiE2dN/{d^Pid^P2) 

^[Pi,P2) {^EidN/d^pi){E2dN/d^p2)' ^^ 

The single- and two-particle cross sections are normalized to the average number 
of particles per event (A^) and the average number of particles in pairs {N{N — \)), 
respectively. Different normalizations of the correlation function are sometimes 
used in the literature |^,^ . Careful consideration of the normalization is required 
when analyzing multi-particle symmetrization effects [0-ll|. 
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2.1.1 Pure quantum statistical correlations 

The most direct connection between the measured two-particle correlations in 
momentum space and the source distribution in coordinate space can be estab- 
lished if the particles are emitted independently ("chaotic source") and propagate 
freely from source to detector. Several approaches to derive this connection |12| 
are worked out in the literature, parametrizing the source as a covariant super- 
position of classical currents p,13-p^ or using a superposition of nonrelativistic 
wave packets |^,|lO|. One finds the simple relations (with the upper (lower) sign 
for bosons (fermions)) 



E — 
C{q,K) 



d X S{x,p) , 



1± 



|/d^x5(x,K)e*9-^| 



jd^xS{x,K + \q) JdSS{y,K 



(9) 
(10) 



where the emission function S{x, K) is an effective single-particle Wigner phase- 
space density of the particles in the source. (Wigner densities are real but not 
necessarily everywhere positive.) For the single-particle spectrum (|9|) this Wigner 
function must be evaluated on-shell, i.e. at p^=Ep={m?+p^Y''^ . The correla- 
tion function (10) was expressed in terms of the relative momentum q=Pi—p2, 



q^=Ei—E2, and average (pair) momentum K={pi+p2)/2, K^=(Ei+E2)/2. As 
the two measured particles are on-shell, p^ 2=-£'l,2=(?7^^+Px 2)^ > ^^^ 4-momenta 
q and K are off-shell. They satisfy the orthogonality relation 



q-K = 0. 



(11) 



That on the rhs of Eq. ( [lO|) one needs the emission function for off-shell momenta 
K may at first seem troublesome. In practice, however, the on-shell approxi- 
mation K^ « EK={'rn?+K'^Y''^ is very accurate in heavy-ion collisions: the 
corrections of order q'^/{8E'j^) are small in the interesting domain of small rel- 
ative momenta q, as a result of the large source sizes and the rest masses of 
the measured particles. A further simplification is achieved by making in the 
denominator of ( p!o[ ) the smoothness approximation |19,2C|, taking the product of 
single-particle spectra at their average momentum K: 



C{q,K) 



1± 



/ d^x S{x, K) e 



iq-x 



Jd^xS{x,K) 



1± 



{e'i-'')iK) 



(12) 



It is exact for exponential single-particle spectra, with corrections proportional 
to their curvature in logarithmic representation. Both approximations can a 
posteriori be corrected for in the correlation radii ("HBT radii", see below), 
using information from the measured single-particle spectra [^. For heavy-ion 
collisions, such corrections are usually negligible |pC|, 21|. 

We call S{x, K) an effective single-particle Wigner density since different deri- 
vations of the relation (^) yield different microscopic interpretations for S. For 
a detailed discussion of this point we refer to |^2|. The differences can become 
conceptually important in sources with high phase-space density |11|. So far, 
in heavy-ion collisions the phase-space densities at freeze-out appear to be low 
enough to neglect them |23,p3]. 
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2.1.2 The invertibility problem 



The mass-sheh constraint ( pAl) ehminates one of the four components of the rel- 
ative momentum q; for example, it can be resolved as 



/3g: 



f3 = K/K" w K/Ek 



(13) 



which gives the energy difference q^ in terms of q and the velocity /3 of the pair. 
With only three independent g-components, the Fourier transform in (|l^ ) cannot 
be inverted, i.e. the space-time structure of S{x,K) cannot be fully recovered 
from the measured correlator: 



C{q,K)-l 



/^e^9-(^-^*)5(x,i^) 



LS{x,K) 



J^e^i-'^Sit,x + pt;K) 



LS{x,K) 



(14) 



Separation of the spatial and temporal structure of the source thus requires ad- 
ditional model assumptions about S{x,K). 

We can connect (14) with (j^) by introducing the normalized relative distance 
distribution which is a folding of the single-particle emission function with itself: 



d{x,K) = I s{X + ^,K)s{X -^,K), s{x',K) 



S{x',K) 



(15) 



L'S{x',K)- 
d is an even function of x. It allows to rewrite |22| the correlator in the form (^): 

C(q, K)-l'^± I d^x cos(g • x) d{x, K) = ± f d^x cos(q • x) Sk{x) . (16) 

In the second equation we used ( [I^ ) and introduced the relative source function 

Sk{x) = I dt d{t, X + j3t; K, Ek) ■ (17) 

In the pair rest frame where (3 = 0, Sk{x) is the time integral of the relative 
distance distribution d, and the time structure of the source is completely inte- 
grated out. On the other hand, Sk{'x) is, for each pair momentum K, fully re- 
constructible from the measured correlator C(q, K) by inverting the last Fourier 



transform in (16). This "imaging method" was recently exploited in |25|. As we 
will see, interesting information about the source dynamics and time structure 
is then hidden in the -fC-dependence of S'x(^); the latter can, however, not be 
unfolded without additional model assumptions about the source. 

2.1.3 Final state interactions and unlike particle correlations 

HBT measurements in high energy physics are mostly performed with charged 
particles. These suffer long-range Coulomb repulsion effects on the way from the 
source to the detector which, even for boson pairs, cause a suppression of the 
measured correlator at q = 0. Moreover, the charged particle pair feels the total 
electric charge of the source from which it is emitted. Final state effects from 
strong interactions play a dominant role in proton-proton correlations |26|, due to 
the existence of a strong s-wave resonance in the two-nucleon channel just above 
threshold. Such final state interaction (FSI) effects are sensitive to the average 
separation of the two particles at their emission points and thus also contain 
relevant information about the source size |4,26, 2^. 
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This has recently led to an increased effort to understand and exploit FSI- 
induced two-particle correlations which also exist between pairs of unlike particles 
[p^32|. The particular interest in such correlations arises from the fact that, 



for particles with unequal masses, they allow under certain circumstances to 
determine the sign of the time difference between the emission points of the 
two particles or the direction of their separation at emission |28-32|; this is not 
possible with correlations between identical particles. In most practical cases, 
however, the FSI- induced correlations are considerably weaker than the Bose- 
Einstein correlations between pairs of identical particles. 



At the level of accuracy of Eqs. (12,16) which use the smoothness approxima- 
tion, the correlator can be easily corrected for 2-body final state interactions by 
replacing e*^'^' with a suitable distorted wave. Instead of ([iq) one thus obtains ||2^ 

(18) 



C{q,K)= j d\SK{r) \^q/2{r) 



A slightly more general result which avoids the smoothness approximation was 



derived in [33|. For simplicity the integral in ( |18D is written in coordinates of the 
pair rest frame where K = Q. ^q/2{f) is an FSI distorted scattering wave for the 
relative motion of the two particles with asymptotic relative momentum q; for 
Coulomb FSI it is given by a confluent hyper geometric function: 



^S/t(^ 



z± 



r(l + r]) e~ 2'^'' 62"-'' F{-ir]; 1; z. 
^{qr ± q ■ r) = ^qr{l it cosO), rj = 



am 



q 



(19) 
(20) 



where a = e^/47r. It describes the propagation of the particle pair from an initial 
separation r in the pair rest frame, at the time when the second particle was 
emitted |33|, under the influence of the mutual FSI. For identical particle pairs it 
must be properly symmetrized: ^q/2 '— > "75(^9/2='= '^-9/2)- For a pointlike source 
Sk{x) = S{x) the correlator (|l^ ) with (^) reduces to the Gamov factor G{r]) 
(to 2G{r]) for identical particles): 

2 _ 27rr] 

For Coulomb FSI it was recently shown |^4| that a very good approximation 
for the Coulomb correction can be taken from measured unlike-sign particle pairs 
in the following form: 



G{ri) = T{l + ir])e 



-\-K-q 



(21) 



cfotAq.K) 



corneas. (Qi ^) Corneas (Qi ^) 



(22) 



G(r?)G(-7?) 

The denominator (which deviates from unity for small q < 8ma) corrects for the 
fact that even for a pointlike source the like-sign and unlike-sign Coulomb corre- 



lations are not exactly each other's inverse. The important observation in |34] is 
that this correction is essentially independent of the source size. 

The effects of the central charge of the remaining flreball on the charged particle 
pair were studied in |27,^. For a static source it was found that at large 
pair momenta the FSI reduces (increases) the apparent size (HBT radius) for 
positively (negatively) charged pairs |2^, ^ whereas for small pair momenta 
the apparent radius increases for both charge states |35]. Expanding sources 
have not yet been studied in this context, nor has this effect been quantitatively 
confirmed by experiment. Also, combining the central interaction with two-body 
FSI remains an unsolved theoretical challenge. 
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2.2 Source Sizes and Particle Emission Times from HBT Correlations 

The two-particle correlation function is usually parametrized by a Gaussian in the 
relative momentum components. We now discuss different Gaussian parametriza- 
tions and establish the relationship of the corresponding width parameters (HBT 
radii) with the space-time structure of the source. 

2.2.1 HBT radii as homogeneity lengths 

The space-time interpretation of the HBT radii is based on a Gaussian approxi- 
mation to the space-time dependence of the emission function S{x, K) [^,36-^. 



Characterizing the effective source of particles of momentum K by its space-time 
variances ("rms widths") 

{x^Xi,){K) = {{x -x)^{x-x)„), (23) 

where (...) denotes the (-ff -dependent) space-time average over the emission func- 
tion defined in (|12|) and x{K) = (x) is the center of the effective source, one ob- 
tains from ( p!^ ) the following generic Gaussian form for the correlator |36,37|,|4H: 



C{q, K) = l± exp [-q^q,{S:f'x''){K)] . (24) 

This involves the smoothness and on-shell approximations discussed in section 
|2.1.l| which permit to write the space-time variances (x^Xiy) as functions of K 
only. Eq. (^) expresses the width of the correlation function in terms of the 
rms widths of the single-particle Wigner density S{x, K). Note that the absolute 
space-time position x{K) of the source center does not appear explicitly and thus 
cannot be measured. 

Instead of the widths of the single-particle function S{x,K) we can also use 
the widths of the relative distance distribution d{x,K) (see (p^)) to characterize 
the correlation function. Starting from (|l^) one finds within the same Gaussian 

approximation C{q,K) = 1 it exp —^q^qi,{x^x'^)d{K) ; here (. . .)d denotes the 
average with the relative distance distribution d. Since d is even, {x^)(i = 0. One 
sees that the rms widths of S and d are related by a factor 2: {x^x")d = 2{x^x^)s- 
This shows that for a Gaussian parametrization of the correlator according to 
(p^), without a factor ^ in the exponent, the width parameters are directly re- 
lated to the rms widths of the single-particle emission function S whereas a similar 
parametrization which includes a factor ^ in the exponent gives as width param- 
eters the rms widths of the relative distance distribution d (or of the relative 
source function 5'j<-(r)). While the latter interpretation may be mathematically 
more accurate, the former is more intuitive and has therefore been preferred in 
the recent literature. 

In either case, the two-particle correlator yields rms widths of the effective 
source of particles with momentum K. In general, these width parameters do 
not characterize the total extension of the collision region. They rather measure 
the size of the system through a filter of wavelength K. In the language intro- 



duced by Sinyukov 1 42 ] this size is the "region of homogeneity" , the region from 
which particle pairs with momentum K are most likely emitted. The space-time 
variances {x^Xu) coincide with total source extensions only in the special case that 
the emission function shows no position-momentum correlations and factorizes, 
Six,K) = f{x)g{K). 
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2.2.2 Gaussian parametrizations and interpretation of HBT radii 



Relating (^4|) to experimental data requires first the elimination of one of the four 



g-components via the mass-shell constraint (11). Depending on the choice of the 
three independent components different Gaussian parametrizations exist. 

A convenient choice of coordinate axes for heavy-ion collisions is the os/-system 
[^,^, with / denoting the longitudinal (or xi) direction along the beam, o the 
outward (or Xq) direction along the transverse pair momentum vector K±, and s 
the third Cartesian direction, the sideward (or Xs) direction. In this system the 
sideward component Ps of the pair velocity (3 in ( p!3| ) vanishes by definition. 

The Cartesian parametrization p4| of the correlator (often referred to, 
historically somewhat incorrectly, as Pratt |13l - Bertsch ^^ parametrization) is 



based on an elimination of g in (24) via dl 



C{q, K) = l± exp 



KiK) 



QiQj 



(25) 



i,j=o,s,l 



The Gaussian width parameters (HBT correlation radii) Rij of the Cartesian 
parametrization are related to the space-time variances of the emission function 
by [|6H|,||] 



Rlj{K) = {{xi- I3ii){xj - (3ji)) , i,j = o, s, /. 



(26) 



These are 6 functions of three variables: the pair rapidity Y, the modulus K± 
and the azimuthal angle $ between the transverse pair momentum K± and the 
impact parameter b. Only these 6 combinations of the 10 independent space-time 
variances {x^x^,) can be measured. 

For azimuthally symmetric collision ensembles the emission function has a re- 
flection symmetry Xs -^ —Xg, eliminating 3 of the 10 space-time variances, and 
the correlator is symmetric under q^ — > —Qs |37]. Then R'^^ = i?^; = 0, and the 
correlator is fully characterized by 4 functions of only two variables K± and Y: 



C{q, K) = l± exp -Riqt - R^ - Rfqf - 2Riiqoqi 



with 



rUk^,y) 

Rf{K±,Y) 






m 



A2\ 



RliK±,Y)-- 
Rli{K±,Y) 



= {{xo- P±i){xi 



(27) 



Pit)). (28) 



paper |^, was recently rediscovered |36 
can be found in Refs. |18,pl, 36-41 



These "HBT radii" mix spatial and temporal information on the source in a non- 
trivial way, and their interpretation depends on the frame in which the relative 
momenta q are specified. Extensive discussions of these parameters (including the 
cross-term R'^^ which originally appeared in the important, but widely neglected 

and then experimentally confirmed p3| , 
18 -p^. The cross-term vanishes in any 
longitudinal reference frame in which the source is symmetric under xi — > — x; |44| 
(e.g. for pion pairs with vanishing rapidity in the center-of-momentum system 
(CMS) of a symmetric collision); in general it does not vanish for pion pairs 
with non-zero CMS rapidity, not even in the longitudinally co-moving system 
(LCMS H) ||. 

For azimuthally symmetric collisions no direction is distinguished for pairs with 
K± = 0. As long as for K± -^ the emission function reduces to an azimuthally 
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symmetric expression (an exception is a certain class of opaque source models 
discussed in section p. 3. 4 ), one has at K_\_ = the identities [^ (x^ — xl) = 
{xqXi) = (ixo) = 0; these imply that R^ — R^ and the cross-term R^^ vanish 
at K± = 0. At non-zero K± these identities for the space-time variances may 
be broken by transverse position-momentum correlations in the source, as e.g. 
generated by transverse collective flow. If the latter are sufficiently weak the 
leading /C^-dependence of the difference 

Rlf, = Rl- Rl = 0i (?) - 2/3^ (x^t) + {xl - ~xl) (29) 

is given by the explicit /3_L-dependence of the first term on the rhs. This yields 
the duration of the particle emission process At = V'(t^) — (t)^ for particles with 
small K^ [p^-|55|. (This is sometimes loosely called the "lifetime" of the effective 
source, but should not to be confused with the total time duration between 
nuclear impact and freeze-out which is not directly measurable.) 

The possibility to extract the emission duration from correlation measurements, 
pointed out by Bertsch and Pratt [^^55|, provided the main motivation for 



the construction of second generation experiments to measure high quality, high 
statistics correlation functions. Subsequent model studies for relativistic heavy- 



ion collisions |52,5q-p9[ where the emission duration is expected to be relatively 
short (of the order of the transverse source extension) showed, however, that the 
extraction of At is somewhat model dependent; the relative smallness of the last 
two terms in (^9|) cannot always be guaranteed, and their implicit iCj^ -dependence 
can mix with the explicit one of the interesting first term. 

The Yano-Koonin-Podgoretskii (YKP) parametrization is an alterna- 
tive Gaussian parametrization of the correlator for azimuthally symmetric col- 
lisions. It uses the mass-shell constraint ( [ill) to express (|J) in terms of q± = 
V^Tql, qi and q^ |^,0,|1|^: 



C{q,K) = l±exp [-Rlql - R^ {qf - (q'f) - {r^ + Rf^ {q ■ Uf] . (30) 

Like in% it has 4 (il'_|_,y)-dependent fit parameters: the three radius parameters 
R±{K), Ru{K), Rq(K), and a velocity parameter U{K) with only a longitudinal 
spatial component: 

C/(K)=7(i^)(l,0,0,r;(K)), 7 = (1 - t;2)~V2. (31) 

The advantage of fitting the form (^) to data is that the extracted YKP radii 
i?_L, R\\ , Rq do not depend on the longitudinal velocity of the measurement frame, 
while the fourth fit parameter v{K) is simply boosted by that velocity. The frame 
in which v{K) = is called the Yano-Koonin (YK) frame; the YKP radii are 
most easily interpreted in terms of coordinates measured in this frame 



Rl{K) = Rl{K) = {xl), (32) 

Rl{K) = {{ii-{(3i/P^)xo)')-{(3i/P±f{il)^{xh, (33) 

RliK) = {(i-Xo/(3^f)-{xl)/l3l^{P), (34) 

where the approximations in the last two lines are equivalent to dropping the last 
two terms in ( |29| ) (see discussion above). To the extent that these hold, the three 
YKP radii thus have a straightforward interpretation as the transverse, longitu- 
dinal and temporal homogeneity lengths in the YK frame. In particular the time 
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structure of the source only enters in Rq. For sources with strong longitudinal 
expansion, like those created in relativistic heavy-ion collisions, it was shown in 
extensive model studies pl,p2,p^59| that the YK velocity v{K) very accurately 



reflects the longitudinal velocity at the center x{K) of the "homogeneity region" 
of particles of momentum K. The YK frame can thus be interpreted as the rest 
frame of the effective source of particles with momentum K, and the YKP radii 
measure the transverse, longitudinal and temporal size of this effective source in 
its own rest frame. 



The parametrizations (27) and ( |30[) use different independent components of q 



but are mathematically equivalent. The YKP parameters can thus be calculated 
from the Cartesian ones and vice versa |41|. The corresponding relations are 

Rl = Rl, (35) 

Rl^ = Rl-R^ = pW[Rl + v^Rf^ , (36) 

Rf = {l-f3f)Rl + ^^l3i-vf[Rl + Rf^, (37) 

Rli = P±{-mii+l''{l3i-v)(Rl + Rf^)) , (38) 

whose inversion reads 



Rl-- 


= B- 


-vC, 


Rq 


= A 


-vC, 


V = 


j-i. -f- u 


2C 


A = 


1 

/3^ 


-^diff) 


B = 


Rf- 




1^, 


-^diff) 



'^ • ■ (38) 



A + B 






C = -—Rli+^Rl^. (40) 



These last definitions hold in an arbitrary longitudinal reference frame. According 
to ( p9| ) V is zero in the frame where C vanishes. However, (p^) also shows that the 
YKP parametrization becomes ill-defined if the argument of the square root turns 



negative. This can indeed happen, in particular for opaque sources [58,61 ; this 
has motivated the introduction of a modified YKP parametrization in 52, ^j, 59 1 
which avoids this problem at the expense of a less intuitive interpretation of the 
modified YKP radii. These remarks show that these relations provide an essential 
check for the internal consistency of the Gaussian fit to the measured correlation 
function and for the physical interpretation of the resulting HBT parameters. 

2.3 Collective Expansion and K -Dependence of the Correlator 

If the particle momenta are correlated with their emission points ( "x-p-corre- 
lations"), the space-time variances in (|2J) depend on the pair momentum K. 
Various mechanisms can lead to such correlations; the most important one for 
heavy-ion collisions is collective expansion of the source. Recently a major effort 
has been launched to extract the collective flow pattern at freeze-out from the K- 
dependence of the HBT parameters. However, thermalized sources may exhibit 
temperature gradients along the freeze-out surface which cause additional x-p- 
correlations. Moreover, pion spectra receive sizeable contributions from the decay 
of unstable resonances some time after freeze-out. These decay pions tend to come 
from a somewhat larger space-time region than the directly emitted ones and, due 
to the decay kinematics, they preferentially populate the low-momentum region. 
Together these two effects also generate x-p-correlations for the emitted pions 
even if the original source did not have them p^ . 
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The separation of these different effects requires extensive model studies some 
of which will be reviewed below. For didactical reasons we will discuss them 
in the context of the YKP parametrization where certain mechanisms can be 
demonstrated most transparently. A translation for the Cartesian fit parameters 
via the cross-check relations ([35|)-(38) is straightforward. Furthermore we will 



show that for sources with strong longitudinal expansion the YK frame (effective 
source rest frame) is usually rather close to the LCMS (in which the pairs have 
vanishing longitudinal momentum), i.e. v{K) ~ in the LCMS. This allows, 
at least qualitatively, for a rather direct extraction of source properties in its 
own rest frame from the Cartesian HBT radius parameters in the LCMS. This is 
important since initially most multi-dimensional HBT analyses were done with 
the Cartesian parametrization in the LCMS, before the YKP parametrization 
became popular. 

We will concentrate on the discussion of azimuthally symmetric sources (central 
collisions) for which extensive knowledge, both theoretical and experimental, has 
been accumulated in the last few years. Many analytical model studies [pi|,|36|-^, 
48,^, 52,57-59, 33,64] are based on the following parametrization of the emission 



function (or slight variations thereof): 



iW^cosh(^7-lO 
b[x,K)- g^^^^ exp 



K-u{x) {r-Tof r^ (??-%) 



T{x) 2(At)2 2R'^{ri) 2{/\r]f 



■ (41) 



Here r'^=x^+y^ , ry=2 ln[(t+2)/(t— z)], and T={t^—z'^Y''^ parametrize the space- 
time coordinates x^ , with d^x = t dr drjr dr d(p. y=2 ln[(-E'i^+i^L)/(£'i<-— -fC/,)] 
and M±={m'^+Kj_y/^ parametrize the pair momentum K. \f2R is the 2-di- 
mensional transverse rms radius of the source ((r^) = 2R'^); R is usually taken 
as a constant, tq is the average freeze-out proper time. At the mean duration 
of particle emission, and Ar/ controls the longitudinal size of the source, L ~ 
TQsinhrj. Note that S{x,K) describes the phase-space distribution at freeze-out 
and not the dynamical evolution of the source from initial conditions; the latter 
is described by hydrodynamical or microscopic kinetic models, discussed below. 
The Boltzmann factor exp[—K-u{x)/T(x)] parametrizes the momentum-space 
structure of the source in terms of a collective, directed component, given by a 
flow velocity field u'^{x), and a randomly fluctuating component, characterized 
by an exponential spectrum with local slope T{x), as suggested by the shape of 
the measured single-particle spectra. Although this parametrization is somewhat 
restrictive because it implies that the random component is locally isotropic, it 
does not require thermalization of the source at freeze-out. But if it turns out 
that all particle species can be described simultaneously by the emission function 



(41), with the same temperature and velocity fields T{x), u{x), this would indeed 
suggest thermalization as the most natural explanation. 
It is convenient to decompose u^{x) in the form 



n^ 



= (cosh rji cosh rjt, sinh rjt e^, sinh r]i cosh r/j 



i^{x) = (cosh ?7; cosh r/t, sinhr/jej,, sinh ry; cosh r/j) , (42) 



with longitudinal and transverse flow rapidities r]i[x) and r]t{x). A simple boost- 
invariant longitudinal flow 7]i{T,7],r) = rj {vi = z/t) is commonly assumed. For 
the transverse flow rapidity profile the simplest choice is 

7]t{T, r/, r) = r]f{T, rj) (^J-^) (43) 
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with a scale parameter r]f. (In our notation, ryj denotes the transverse collective 
flow rapidity, whereas f3± is the transverse velocity of the particle pair.) In most 
studies r]f was (like R) set constant such that r]t{r) was a function of r only. This 
cannot reproduce the observed rapidity dependence of {p±) and of the inverse 
slopes of the m_L-spectra [^|6^. In [^^ it was shown that an ry-dependence 
of R, with R{rj) shrinking in the backward and forward rapidity regions keeping 
the slope rjf/R in ( [43| ) fixed, is sufficient to fully repair this deficiency. Here we 
discuss only the simpler case of constant R, r]f . 

With these assumptions the exponent of the Boltzmann factor in (^l|) becomes 



K • u{x) 



= M± cosh(y — rj) cosh rjt{r) — K±-er sinh?7((r) . 



(44) 



For vanishing transverse fiow {r]f = 0) the source depends only on M_|_ and 
remains azimuthally symmetric for all K±. Since in the absence of transverse 
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Figure 2: The YKP radii R±, R\i, and Rq (top to bottom) as functions of M± 
for pairs at Ycm = 0. Left: no transverse fiow. Right: rjf = 0.6. Solid (dashed) 
lines are for pions (kaons). Note the breaking of M_L-scaling by transverse fiow. 
Source parameters: T = 140 MeV, Ary = 1.2, R = 3 fm, tq = 3 fm/c. At = 1 
fm/c. (Figure taken from ||49||.) 



flow the /3-dependent terms in (^3|) and (|3J) vanish and the source itself depends 
only on M±, all three YKP radius parameters then show perfect M±- scaling: 
plotted as functions of M_\_ , they coincide for pion and kaon pairs (see Figure S, 
left column). This remains true if T{x) varies with x; temperature gradients in 
the source do not destroy the M_L-scaling. 
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For r]f ^ (right column) this M_|_-scaling is broken by two effects: (1) The 
thermal exponent (|43) receives an additional contribution proportional to K±. 



(2) The terms which were neglected in the approximations (33,^) are non-zero, 
and they also depend on /3± = K±/Ek- Both induce an explicit rest mass 
dependence and destroy the Mj^-scaling of the YKP size parameters. 

2.3.1 Longitudinal flow: Yano-Koonin rapidity and M_|_ -dependence of Rn 

At each point in an expanding source the local velocity distribution is centered 
around the average fluid velocity u{x). Thus two fluid elements moving rapidly 
relative to each other are unlikely to contribute particles with small relative mo- 
menta. Only source regions moving with velocities close to that of the observed 
particle pair contribute to the correlation function. How close is controlled by the 
width of the random component in the momentum distribution: the larger the 
local "thermal smearing" , the more the differences in the fluid velocities can be 
balanced out, and the larger the "regions of homogeneity" in the source become. 
Longitudinal expansion is most clearly reflected in the behaviour of the Yano- 
Koonin (YK) rapidity Yy-^=^ln[{l+v)/{l—v)]. Figure ^ shows (for pion pairs) 
its dependence on the pair momentum K. Transverse flow is seen to have a 
negligible influence on the YK rapidity. On the other hand, the linear dependence 
of Yyj^ on the pair rapidity Y (Figure ya) is a direct reflection of the longitudinal 



expansion flow |41|; for a non-expanding source Yy^ would be independent of Y. 
The correlation between the velocities of the pair (Y) and of the emission region 
(Yyk) strengthens as the thermal smearing decreases. For the Boltzmann form 
(pl|) the latter is controlled by T/M± , and correspondingly Y^k (-^-L ) ^) ^ ^ 
as K± — > oo. For small K±, thermal smearing weakens the correlation, and 
the effective source moves somewhat more slowly than the observed pairs, whose 
longitudinal velocities have an additional thermal component. 

If the ratio of the longitudinal source velocity gradient to the thermal smearing 
factor, defined in (^) below, is large, Rh becomes small and the longitudinal 
rapidity Y^k o^ ^^^ effective source becomes equal to that of the emitted pairs, Y. 
This can be true even for slow longitudinal expansion as long as it is strong enough 
compared to the thermal smearing. Consequently, observation of a behaviour 
like the one shown in Figure ^ demonstrates strong, but not necessarily boost- 
invariant longitudinal flow. 

Longitudinal expansion is also reflected in the M_L-dependence of Rn (second 
row of Figure |2|). Comparison of the left with the right diagram shows only 
minor effects from transverse expansion; longitudinal and transverse dynamics 
are thus cleanly separated. A qualitative understanding of the M_|_-dependence 
is provided by the following expression, valid for pairs with Y = 0, which can be 



derived by evaluating (|33D via saddle-point integration |21,38, 3£, S7|: 



2 



Ru^Ll = ^^2w (45) 



Eq. (03) shows explicitly the competition between the longitudinal velocity gra- 
dient d-ui and the thermal smearing factor T/M±. For strong longitudinal ex- 
pansion (large velocity gradient and/or weak thermal smearing) and/or large 
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geometric longitudinal extension tqAij of the source the second term in the de- 
nominator can be neglected, and i?|| drops steeply as l/^/M± |67|. Note that 
quantitative corrections to (^) are not always small [pC[|. 

We emphasize that only the first equation in (|46|) is general. The appearance 



of the parameter tq in the second equation is due to the choice of a Bjorken 
profile for the longitudinal flow for which the longitudinal velocity gradient is 
given by the total proper time tq between impact and freeze-out. This is not true 
in general; the interpretation of the length Ru in terms of the total expansion 
time is therefore a highly model-dependent procedure which should be avoided. 
As a matter of principle, the absolute temporal position of the freeze-out point 



is not measurable, see section 2.2.1 
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Figure 3: (a) The Yano-Koonin rapidity for pion pairs, as a function of the pair 
cm. rapidity y, for various values of K^ and two values for the transverse flow 
?7/. (b) The same, but plotted against K^ for various values of Y and t]j. Source 
parameters as in Fig. y. (Figure taken from [41|.) 



Strong longitudinal x-p-correlations also occur in string fragmentation. In fact, 
in the Schwinger model of string breaking |£^] the quark pairs created from 
the chromoelectric field are assumed to have longitudinal momentum Y = r\^ 
without thermal fluctuations. Thus a similar linear rise of the YK-rapidity with 
the pair rapidity and a strong decrease of R\\{Mx_) would also be expected in 
jet fragmentation (with the x^-axis oriented along the jet axis). It would be 
interesting to confirm this prediction p^ in e^e~ or pp collisions. 



2.3.2 Transverse flow: M_L-dependence of R^ 

Just as longitudinal expansion affects R\\ , transverse flow causes an M_L-depen- 
dence of R±. This is seen in the first row of diagrams in Figure §, which also 
shows that longitudinal flow does not contribute to this feature: for ijf = the 
transverse radius does not depend on M±, in spite of strong longitudinal expan- 
sion of the source. A qualitative understanding of this behaviour is given by 
the analogue of (|45|), again obtained by evaluating (32) via saddle point integra- 



tion |2l|,|33,| 



Ri 



Rl^ 



i?2 



r:' 



1 + {R/R^y 1 + VJiM^/T) 



inflow (Mx) 



dr]t{r)/dr y M± 




(47) 



(48) 
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Once again there is competition between flow velocity gradients in the source, this 
time in the transverse direction, which tend to reduce the homogeneity regions, 
and thermal smearing by the factor T/M_\_ resulting from the random compo- 
nent in the momentum distribution, enlarging the regions of homogeneity. The 
left equations in (4^,4^) are generic while the right ones apply to the specific 



transverse flow profile (p 

Transverse fiow must be built up from zero during the collision while longi- 
tudinal x-p-correlations may contain a sizeable primordial component from in- 
complete stopping of the two nuclei and/or the particle production process (e.g. 
string fragmentation, see above). One thus expects generically weaker transverse 
than longitudinal fiow effects at freeze-out. Correspondingly, in realistic simula- 
tions (e.g. [^, Qg,|5^) the longitudinal homogeneity length i?|| turns out to be 
dominated by the expansion (i.e. by Lflow) while in the transverse direction the 
geometric size R dominates at small M_|_, with fiow effects taking over only at 
larger values of M_|_. Correspondingly, R±{M_i) decreases more slowly at small 
M^ thani2||(M^) [||,|7|- 

2.3.3 The emission duration 

Saddle-point integration of (34) with the source (^Tj) yields, with L^ from (p5|). 



i?2«(At,)2^(Ar)2 + 2K/r2 + L2-ro . (49) 



2 



The M_L -dependence of L^, thus induces an M_L-dependence of the temporal YKP 
parameter. Eq. (^) refiects the proper time freeze-out assumed in the model 
([4l|): particles emitted at different points z are also emitted at different global 
times t, and the total temporal width of the effective source is thus given by the 
Gaussian width At plus the additional variation along the proper time hyperbola. 



integrated over the longitudinal homogeneity region L* [40|. 

Although (^9|) suggests that the proper emission duration At can be measured 
via Rq in the limit M_|_ -^ cxd (where L^, — > 0), this is has not been true in systems 
studied to date. For transversely expanding sources Ro receives additional con- 
tributions from the /3_L-dependent terms in (|34|), in particular at large M_l (see 
Figure S). The most important correction is due to the term (x^ — x^)//?^ which 



can have either sign and usually grows with M±_ |^,5£, 61|. The extraction of 
the emission duration must thus be considered the most model-dependent aspect 
of the HBT analysis. 

2.3.4 Temperature gradients and opacity effects 

A different source of transverse x-p-correlations which can compete with trans- 
verse fiow in generating an M_L-dependence of R±_ are transverse temperature 
gradients in the emission function. Since particle densities and mean free paths 
(which control the freeze-out process |69|,|f0|) depend very strongly on temper- 
ature, one would a priori not expect strong temperature variations across the 
freeze-out surface [^. While transverse temperature gradients and transverse 
flow affect i?_L similarly (except that the latter weakly breaks the M_|_-scaling), 
they have quite different effects on the single particle spectra 1 3£ , 5C ] : transverse 



temperature gradients strongly reduce the flattening effect of transverse flow on 
the r?T,_|_-spectra which is needed to reproduce the data [p0|,p2|. Thus constrained 
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by single-particle data, their phenomenological usefulness is limited. Temporal 
temperature gradients only reduce the emission duration, but do not affect the 
transverse i?-parameter p9|,pC|]. 

One possible source feature that parametrization ( p|) cannot describe is "opac- 
ity", i.e. surface dominated emission. Heiselberg and Vischer generated opaque 
sources by multiplying ( p| ) (or a similar source with bulk freeze-out) with an 
exponential absorption factor |6^ (see also 



Sopiiqne{x,p) = 5(x,p)exp -a/S/tT (/cff/Amfp) 




(50) 



^eff = les{r,(l)) = e "^ e ^^ dx' . (51) 



Xo 



The ratio X^afp/R controls the degree of opacity of the source; as Xmfp/R -^ 0, 
the source becomes an infinitely thin radiating shell. The parametrization ( |50|) 
together with (|4l| ) leads to sources with negative (x^ — x'^) for all values of K± 
(including the limit K± -^ 0). According to (|29| ) and (M) this leads to negative 
values of R%g and Rq {Rq even diverges as K± — > [p8[); the data (see below) 
are consistent with vanishing or positive -R^jg at small K±. 

A non- vanishing difference (x^— x^) in the limit K± — > violates the postulated 
azimuthal symmetry of the source (see discussion before Eq. (|29|)). It is easy to 
see that short-lived sources can never be opaque for particles with K± -^ 0: the 
source shrinks to zero before such particles can be reabsorbed. The particular 
behaviour excluded in JS^ is thus anyhow rather unphysical. At larger K±, on 
the other hand, the "opacity signal" (x^— x^) < (leading, if strong enough, 
to i?o < R^ [^) can be "faked" by other mechanisms: Tomasik found \5%, 59| 



that expanding sources with a box-like transverse density profile generate exactly 
such a signature. At the moment it is thus unclear how to uniquely distinguish 
"opaque" from "transparent" sources. 

2.4 Non-Gaussian Features of the Correlator and q-Moments 

The (/f-dependent) HBT radii provide a full characterization of the two-particle 
correlation function only if it is a Gaussian in q or, equivalently, if the effective 
source ^(x, K) is a Gaussian in x. However, in many physical situations the 
source is not characterized by just one, but by several distinct length scales. In 



this case the Gaussian approximation of section 2.2 breaks down. 



2.4.1 Resonance decays 

The most important physical processes leading to a non-Gaussian shape of the 
correlator are resonance decays [18, 32, 73~|77]. Especially longlived resonances 



which decay into pions cause a long-range exponential tail in the pion emission 
function which distorts the two-particle correlation function at small relative 



momentum q (see example in Figure 4 below). According to |77,75] the resonances 
can be classified into three classes: 

1. Short-lived resonances (F > 30 MeV) which (especially if heavy) decay 
very close to their production point. Their most important effect is to 
add a contribution proportional to their lifetime to the emission dura- 
tion [^, 77|, thereby affecting i?^ and R^ in the Cartesian and iJg in the 
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YKP parametrization, but not the transverse radius R±. They do not spoil 
the Gaussian parametrization. 

2. Long-lived resonances (T ^ 1 MeV), mostly the r],r]',K'^ and hyperons. 
These resonances travel far outside the original source before decaying. The 
resulting wide tail in the emission function contributes to the correlator only 
at very small relative momenta. This region is experimentally inaccessible 
due to finite two-track and momentum resolution, and the contribution from 



this "halo" [76| to the correlator is thus missed in the experiment. The result 
is an apparent decrease of the correlation strength (i.e. the intercept A). In 
the measurable g-range the shape of the correlator is not affected. 

3. The uj meson (F = 8.4 MeV) is not sufficiently longlived to escape detection 
in the correlator. Its lifetime is, however, long enough to create a measurable 
exponential tail in the pion emission function which distorts the shape of 
the correlator, giving it extra weight at small q and destroying its Gaussian 
form. 

In practice the pions from short-lived resonances can thus be simply added 
to the directly emitted ones into an emission function for the "core" ||73| , ^ . 
The "halo" from long-lived resonances is accounted for by a reduced intercept 
parameter 

\{K)=U-Y,fr{K)\ , (52) 

where the sum goes over all longlived resonances and fr{K) is the fraction of pions 
with momentum K stemming from resonance r. A correspondingly modified 
Cartesian parametrization for the correlator reads 

C{q, K) = l± X{K) exp I - ^ RU^)mj J • (53) 

Pions from to decays must, however, be considered explicitly and, if sufficiently 
abundant, the resulting correlator is no longer well described by the ansatz (]5^). 
In heavy-ion collisions the resonance fractions fr are unknown since most reso- 
nances cannot be reconstructed in the high-multiplicity environment. Thus X{K) 
in (^3[) is an additional fit parameter. Its value is very sensitive to non-Gaussian 
distortions in the correlator, and so are the HBT radii extracted from a fit to the 



function (|53|). In theoretical studies |77] it was found that differences of more 



than 1 fm in the fitted HBT radii can occur if the fit is performed with A fixed 
to its theoretical value ( p2|) or if, as done in experiment, A is fitted together 
with the radii. In the latter case resonance contributions (including the uj) af- 
fect the fitted radii much less than in the former. This difference in procedure 
may largely explain the consistently larger resonance contributions to the HBT 
radii found by Schlei et al. |73,75|, compared to the much weaker effects reported 



in |77|. Whereas Schlei et al. |73, 7^ find that resonances, whose decay pions con- 
tribute only to the region of small K±, add considerably to the M_|_-dependence 
of Rg = R± and thus contaminate the transverse flow signature, practically no 
such effect was found in [77|. 
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2.4.2 Q'-moments 

In view of these systematic uncertainties one may ask for a more quantitative 
characterization of correlation functions whose shape deviates from a Gaussian. 
This can be achieved via the so-cahed g-moments of the correlator [^]. In this 
approach the matrix of Cartesian HBT parameters TZ = (RfAK)) (i,j = o,s,l) 
and the correlation strength X{K) are calculated from the following integrals: 

A(ii:) = yd^t^^KJ7^ /"d3g[C(q,ii:)-l]. (55) 

Similar expressions exist for the YKP parameters [Q. For a Gaussian corre- 
lator this gives the same HBT parameters as a Gaussian fit; for non-Gaussian 



correlators the HBT radius parameters and intercept are defined by (54,55). 

Deviations of the correlator from a Gaussian shape are then quantified by 
higher order g-moments. Formally they can be obtained as derivatives at the 
origin of the relative source function SKi^) which acts as generating function ||77|| : 

Practical applications of this method are limited by severe statistical requirements 
for the measured correlator. So far they have been restricted to uni-directional 
moments along one of the three g-axes. The leading deviation of C{q, K) from a 



Gaussian shape is then given by the kurtosis [77| 

Ml 



MK) = ,^f;%^{K)-l, i = o,s,l. (57) 



In 1 77 1 the influence of transverse flow and resonance decays on the transverse 
HBT radius Rg = R± on the kurtosis A^ {K± ) was studied for sources of the type 
(PI). It was found that decay pions give a positive contribution to the kurtosis 
which disappears at large K± together with the resonance fractions fr{K). Trans- 
verse flow, on the other hand, leads to a vanishing or very small negative kurtosis 
which tends to become larger with K±. The sign and iC_L -dependence of the 
kurtosis thus provide a possibility to check whether a measured -fC^-dependence 
of Rg is really due to transverse flow or "faked" by resonance decays fT^] • 

2.5 The Average Freeze-out Phase-Space Density 



Bertsch ||79[] pointed out that by combining measurements of single-particle mo- 
mentum spectra and two-particle correlations one can determine the spatially 
averaged phase-space density at freeze-out and thereby test local thermal equi- 
librium in the pion source created in high energy nuclear collisions: 

_ J^f{x,p)pf^d^a^{x) 
^^^^P^- J^f{x,p)p^d^a,{x)- ^^^> 

Here d^a{x) is the normal vector on a space-like space-time hypersurface S(x). 
According to Liouville's theorem, S is arbitrary as long as its time arguments are 
later than the time tf{x) at which the last pion passing the surface at point x 
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was produced. If the measured single-particle p^-spectrum is parametrized by an 
exponential with inverse slope parameter TgQ^y) and the two-particle correlation 
function by the Gaussian (p^), one finds [p2|,^[79| 



{f){K±,Y) 



^X{K^,Y){dn/dY) (2^rjj(y))-i g-^^/^-ff^^) 



^-3/2 Ep Rs{K^,Y)jRl{K^,Y)Rf{K^,Y) - R%{K^,Y) 



■ (59) 



The numerator (where dn/dY denotes the multiplicity density of a single charge 
state) gives the momentum-space density at freeze-out while the denominator re- 
flects the space-time structure of the source and can be interpreted as its covariant 
homogeneity volume for particles of momentum K. The factor vA ensures [^] 
that only the contributions of pions from the decays of short-lived resonances, 
which happen close to the primary production points, are included in the average 
phase-space density (see section 2.4.1 ). 

2.6 The Usefulness of 3- Particle Correlations 

Two-particle correlations are insensitive to the phase of the two-particle exchange 
amplitude: writing the latter for two particles with momenta pjj as 



d'^xSix, \{p, +p,)) e^(P-P.)- = p,^ e'<t 



(60) 



the phase (/>jj is seen to drop out from the correlator (12). This is no longer true 
for higher-order multiparticle correlations. For example, for a completely chaotic 
source the true 3-particle correlator, with all two-particle correlation contribu- 
tions R2{i,j) = C{pi,pj) — 1 removed, 

R3{Pl,P2,P3) = C3ipi,P2,P3) - i?2(l,2) -i?2(2,3) - i?2(3, 1) - 1, (61) 

and properly normalized, 

R3{Pl,P2,P3) /«„x 

^^^P^^P^^P^^ = Vi^,(l, 2) i^2(2, 3)^2(3,1) ' ^''^ 

gives the sum of phases of the three two-particle exchange amplitudes |8^ : 

r3{Pi,P2,P3) = 2 cos (012 + 023 + 03i) = 2cos$. (63) 

Expanding the two-particle exchange amplitude ( |60D for small relative momenta 
one finds |80| 



$ 



2 '^12 123 



d{x^)3 d{Xu)3 



dK^ 



^K^' 



1 

24 
1 
2 



[<ll2 <f\2 Q23 + Q23 923 Qu] 



d'^{x 



m/3 , d'^{Xu)3 



+ 



+ 



d\xx), 



(ll2 923 (912 + 923)^ {Xf^XiyXxIS + 0{q^) . 



(64) 



Here qij = Pi—pj (with qi2 + q23 + '731 = 0) are the two-particle relative momenta 
and K = (pi + P2 + Ps)/^ is the average momentum of the particle triplet; the 
averages (. . .)3 are calculated with the emission function S{x, K). Eq. (|6^) shows 
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that $ depends on the odd space-time variances (x^) etc. of the emission function 
and on the derivatives with respect to K of the point of highest emissivity (3;) 3. 
These reflect the asymmetries of the source around its center. In the Gaussian 



approximation of section 2.2.1 they vanish. 

These considerations show that the true three-particle correlator contains ad- 
ditional information which is not accessible via two-particle correlations. In prac- 
tice, however, it is difficult to extract this information. The leading contribution 
to $ is of second order in the relative momenta qij , and in many reasonable models 



it even vanishes |81|. Therefore new information typically enters ?'3(pi,P2,P3) at 
sixth order in q. The measurement of the phase $ is thus very sensitive to an 
accurate removal of all leading g^-dependences by a proper determination and 
normalization of the two-particle correlator. 

On the other hand, it was pointed out that the intercept of the normalized true 
three-particle correlation parameter r^ may provide a good test for the chaoticity 
of the source. Writing the emission function for a partially coherent source as 
S = Scha. + "^coh and denoting the chaotic fraction of the single-particle spectrum 
at momentum p as e(p), the intercept of r^ is given by [p0||82|| 



A3(K) ^ r,iK, k, k) = ^^<^) ^2-e{k))l^ - ^^^^ 

This relation is useful since, contrary to the two-particle correlation strength A, 
the intercept ( |65|) of the normalized three-particle correlator is not affected by 



decay contributions from long-lived resonances which cancel in the ratio (62) 



Complete small-g expansions of R2 and R3 which generalize the Gaussian 



parametrization (24) to the case of partially coherent sources and to three-particle 
correlations, improving on earlier results in [S2, 33|, can be found in [^0[. Within 
a multidimensional analysis of 2- and 3-pion correlations they allow separate de- 
termination of the sizes of the homogeneity regions of the chaotic and coherent 
source components as well as the distance between their centers. 

3 TWO-PARTICLE CORRELATIONS FROM DYNAMICAL 
MODELS 

Interpretation of correlation functions measured in heavy-ion collisions requires 
understanding the true relationship of the parameters extracted from fitting the 
data and the actual single-particle distributions at freeze-out. The level to which 
this works in practice can be established by using an event generator to model 
the collision dynamics, particle production and hadronic rescattering, and then 
constructing a two-particle correlation function. These functions can be fit in the 
same way as experimental correlation functions and the fit parameters compared 



to the single-particle freeze-out distribution in the event generator [78,84|. 

The event generator correlation functions are constructed from the positions 
and momenta representing the single particle emission distribution at the time 
of the last strong interaction (i.e. at freeze-out). The subsequent calculation 
of correlation functions uses particle pairs drawn randomly from this list and 
constructs a two-particle symmetrized wave function |22,5^, ^, ^]. Coulomb 
wave functions for the particles are used. As for experimental data, a Coulomb 
correction is applied to the correlation function before fitting [[7^]. 

As the event generator yields a correlation function while simultaneously know- 
ing the space-time distribution and history of the same particles, discrepancies 
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between the fit parameters and freeze-out distributions may be resolved. Fur- 
thermore, the generated particles may be subjected to experimental acceptance 
cuts and treated like real data. This allows evaluation of the effects of experi- 
mental acceptances and analysis techniques. Since a significant number of the 
observed hadrons arise from decay of long-lived resonances, the event generators 
can also quantify their effects on the correlation functions. Such studies were 
performed using the RQMD [78, g6-|89|] and ARC event generators [84, 90| as well 
as hydrodynamical simulations [73,|75|, 91-p4[. 



3.1 RQMD 



Many experiments use the RQMD event generator [36[, as it satisfactorily repro- 
duces single-particle distributions. RQMD simulates the microscopic phase-space 
evolution, using resonance and string excitation as primary processes, followed 
by fragmentation, decays and subsequent hadronic collisions. Many features in 
p-nucleus and nucleus-nucleus collisions which can be related to secondary scat- 
tering are well described by RQMD [p5|,pH[. The numerous secondary collisions 



result in considerable transverse flow of RQMD events before freeze-out [56,78,97[. 



3.1.1 Collective expansion 

Figure ^ shows a comparison of RQMD freezeout positions and correlation func- 
tions for S+Pb collisions at 200 GeV/nucleon |7^. The top half shows the pion 
(solid lines) and kaon (dotted lines) position distributions at freeze-out along and 
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Figure 4: Comparison of position distributions at freeze-out and calculated corre- 
lation functions, both from RQMD events. The lines in the bottom figures show 
1 + [p(g)[^ where p{q) is the Fourier transform of the corresponding distribution 
in the upper figures. The notation on the axes 9par ^ 'Ztside ' £xput corresponds to 
qi,qs,qo used elsewhere in this review. (Figure taken from flSj.) 



transverse to the beam direction as well as in time, integrated over their momenta. 
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The kaon distributions are narrower than those of the pions. The lower section of 
the figure shows correlation functions calculated from the RQMD events, plotted 
as functions of the Cartesian variables. The points show the calculated correlation 
functions using the NA35 experimental acceptance, while the solid and dashed 
lines indicate Fourier transforms of the relevant components of the pion and kaon 
freeze-out distributions in the top row of figures. If the particle positions and 
momenta were uncorrelated, HBT interferometry should reproduce the full size 
of the freeze-out distribution, and in the lower panels the lines should agree with 
the corresponding points. Instead, the correlation functions calculated from the 
RQMD phase-space distribution are much wider in qi than expected from the 
Fourier transform, indicating a smaller effective source. This reflects longitudinal 
position-momentum correlations arising from a strong longitudinal expansion of 
the source. As will be seen later, this prediction by RQMD and hydrodynamical 
models is confirmed by the data. An analogous effect is seen in the sideward 
direction where especially the calculated kaon correlation function is wider than 
expected from the Fourier transform of their momentum-integrated transverse 
spatial distribution. This can be traced back to transverse position-momentum 
correlations in the RQMD freeze-out distribution, induced by collective trans- 
verse expansion. As discussed above, these correlations increase with increasing 
transverse mass of the particles and are thus more strongly reflected in the kaon 
correlation function. The resulting decrease of the transverse effective source size 
extracted from kaon correlations is again confirmed by the data. 

Event generators like RQMD can be used to study the infiuence of experimen- 
tal acceptance cuts on the ir_L-dependence of the HBT fit parameters induced 
by collective flow. RQMD was shown to reproduce the -fC_L -dependence of the 
correlation functions both at SPS 87-|89[ and AGS |98, 9^ energies. It is note- 



worthy that by analyzing the same set of RQMD events with the NA35 and NA44 
acceptances, good agreement with both data sets was obtained even though the 
iir_L-dependence of Rs appears to be somewhat stronger in NA44 [^]. 

The shrinking of the effective source size as a result of collective flow is il- 
lustrated in Figure ^. The collective flow velocity has the effect of "focussing" 
particles arising from nearby regions of the source. As the correlation selects 
particles of small momentum difference, it is sensitive to this focussing. Figure g 
shows the freeze-out distributions of all pions (open histograms) and of those pi- 
ons accepted in two NA44 spectrometer settings (44mr is the low K± setting and 
131mr high K±). Higher K± pairs have a narrower distribution at freeze-out due 
to the larger effect of the flow velocity. Note that this is the focussing from flow 
and not an experimental acceptance effect. The correlation function only "sees" 
particles which are flowing in the same direction; particles from the far side of 
the source flow away in an expanding source. This is why the characteristic drop 
of apparent source size with increasing K± is observed by both NA35 and NA44 
even though one experiment has large acceptance and one narrow. 

3.1.2 Emission duration 

The emission duration is usually calculated from the flt parameters via Eq. (| 
using the leading term only. Comparison of this emission time estimate with the 
actual width of the freeze-out time distribution from RQMD events determines 
the validity of neglecting the position-momentum correlation effects on the space- 
time variances. Fields et al. |56| found that only for K± below ~ 100 MeV/c 
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Figure 5: RQMD freeze-out distributions for pions. The open histograms give 
all pions from RQMD, the hatched ones only those accepted by NA44 at various 
spectrometer settings, x, y, z indicate the out-, side- and longitudinal directions 
in the collision center of mass system. (Figure taken from 



does this method yield the actual emission duration. For higher K± particles, 



the true emission duration (3-6 fm/c for S+Pb |56] and 7-8.5 fm/c for Pb+Pb 
p9| ) is shorter due to flow- induced position- momentum correlations, but the 
value derived from the difference between R^ and R^ underrepresents even this. 
Consequently, extraction of the emission duration from experimental correlation 
functions should only be attempted at low K±. 

Of course, the total lifetime of the source between impact and freeze-out is 
considerably larger. 



3.1.3 Resonance decays 

In RQMD as many as half of the low K± pions in heavy-ion collisions arise from 
the decay of long-lived resonances ||78|. Very long-lived resonances produce a 
long tail in the pion freeze-out position distribution, visible in Figure ^. This 
corresponds to a component of the correlation function too narrow to measure 
and reduces the correlation strength A |^]. In S+Pb collisions the fraction of 
pions from uj,r],rj' resonances is 30% at low p± and falls to 5 % at p±= 800 
MeV/c [^]. These pions cause a departure from a Gaussian source shape, which 
impacts the quality of fits with a Gaussian parametrization. 

Resonance decay contributions to kaons are smaller, which makes the interpre- 
tation of their correlation functions cleaner. Due to their rest mass they always 
have a sizeable M±, and so their correlation functions are more sensitive to flow. 
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3.2 Hydrodynamical Models 

A more macroscopic approach to describe the dynamical evolution of the reac- 
tion zone treats it as a locally thermalized ideal fluid and solves the relativistic 
hydrodynamical equations. Initial conditions are usually set after an initial pre- 
equilibrium stage and suitably parametrized |7^, 75|, pl| , |100| ]. Freeze-out is usually 



enforced at a fixed energy density or temperature. After the end of the simu- 
lation those fluid cells which satisfy this freeze-out criterium are identified, and 
the local energy and baryon density in these cells are converted into temperature 
and chemical potentials using the equation of state of an ideal resonance gas in 
thermal and chemical equilibrium. Each such cell thus emits a thermal hadron 
spectrum boosted by the local fluid velocity. This determines the emission func- 
tion S{x,p) of the model from which spectra and correlation functions can be 



calculated |7|,|7|,||,|3| 



Measured spectra and correlation functions put constraints on the output of 
such simulations which can be used to identify allowed combinations of initial 
conditions and equations of state 1 94, 100(1 . In this way certain classes of evolution 



scenarios can be eliminated while successful combinations can be used to predict 
other observables for further discrimination or hadronic one- and two-particle 



spectra at future colliders |101|. 



4 TWO-PARTICLE CORRELATIONS IN HEAVY-ION EXPE- 
RIMENTS 

4-1 General Remarks and Short Overview of the Experiments 

4.1.1 Construction of the correlation function 

Experimentally, correlation functions are constructed by counting events with 
boson pairs of given pair and relative momenta and dividing by a properly nor- 



malized [|,102| "background" sample with no enhancement: 

CT'"^'-iq,K)=Aiq,K)/B{q,K). (66) 

Typically, B(q,K) is generated by creating artificial pairs by combining single 
tracks from different events, or by using unlike-sign pions. Generally, analyses 
bin both the data and the background in the chosen variables. The correlation 
functions are then corrected for Coulomb interactions, experimental resolution 
and two-particle acceptance (generally by Monte Carlo techniques), residual two- 
particle effects on the single-particle spectrum constructed from mixed pairs, 
particle misidentification (if contamination is significant); see [103-lOf:] for de- 
tails. 

4.1.2 Fit procedures 

Correlation functions are customarily fit by a Gaussian in q which assumes a 
Gaussian-distributed source. This is a simplification, since the source may well 
have a more complex shape. However, the Gaussian assumption provides a rea- 
sonable representation of the data |105 , 107] and continues to be used. One should 



keep in mind, however, that the measured correlation functions are usually not 
perfect Gaussians and that this leads to systematic uncertainties of the order of 
10-20% in the extracted fit parameters. 
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High statistics data samples are now analyzed using multi-dimensional fits 
|13,36,37, 41, 57|. Unlike many lower-dimensional parametrizations these do not 



require the unrealistic assumption of a spherical source and are more sensitive 
to the collision dynamics [^. Some analyses do, however, use lower-dimensional 
parameterizations due to statistical or acceptance limitations. These carry the 
danger of producing misleading results, but they can be useful if they are based 
on a complete parametrization and make proper projections [p^. 

4.1.3 Coulomb corrections in experiment 

Initially, experimental correlation functions were corrected for 2-body Coulomb 
interactions using the Gamov factor. However, as measurements of heavier sys- 
tems were made, the point source approximation became increasingly inappro- 
priate. Several different techniques have been used for improved corrections. 
The first experimental improvement was achieved by integrating expression 



( [Tq ) with Coulomb wave functions using a technique developed by Pratt |105]. 
For the relative source function Si<-('') one usually takes a spherically symmetric 
Gaussian in the pair rest frame with a size parameter which is iterated. In S+Pb 
collisions, NA44 found that the improved Coulomb correction resulted in a 5- 
10% increase in the fitted radius parameters for pions and kaons [gS, 103, 104]. In 



Pb-|-Pb collisions, the difference is 8-12%, with A decreasing by 3-6% [p9 |. 

The Coulomb correction can be investigated experimentally by measuring cor- 
relation functions of oppositely charged particle pairs, where quantum mechanical 
symmetrization effects are absent |109(| . The measured Coulomb attraction re- 



flects the (non-zero) source size |2^, and can be used to parametrize a Coulomb 



correction for same sign pairs. This technique is used by NA49 107 |. 



4.1.4 Short description of the experiments 

We now briefly describe the experiments whose data will be discussed in the fol- 
lowing. E802/859/866 at the AGS is a wide acceptance magnetic spectrometer 
experiment, tracking particles with drift chambers and identifying them via time- 
of-flight. E814/877 at the AGS is a multipurpose experiment which includes a 
hadron spectrometer covering forward angles. Particle identification is achieved 
using time-of-flight with a very long flight path. The angular range of the accep- 
tance limits the identified hadrons to rather small p±. E895 is a time projection 
chamber (TPC) at the AGS, which identifies hadrons through their energy loss 
in the gas-filled detector volume. 

At the CERN SPS, NA44 is a second generation experiment which measures 
single and two-particle distributions at midrapidity. It is characterized by excel- 
lent particle identification, with contaminations at the 1% level. As a focussing 
spectrometer, its acceptance for particle pairs with small momentum difference 
is optimized, allowing for high statistics in the region of the Bose-Einstein cor- 
relation signal. NA35 is a streamer chamber and TPC experiment at CERN, 
with a large acceptance for pions from 200 A GeV S-nucleus collisions. NA49, 
the successor to NA35, uses four TPC's and two time-of-fiight walls to track and 
identify particles in Pb-|-Pb collisions. For the correlation functions NA35 and 
NA49 use charged tracks without particle identification; they are fit as a function 
of pair rapidity and transverse momentum. 

Different experiments analyze correlations in different reference frames. The 
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NA44 Collaboration |88, 104 1 use the LCMS in which the longitudinal pair mo- 



mentum vanishes. This frame couples the lifetime information solely to Qq and 
ensures that the source velocity in the analysis frame is usually small [^. The 



NA35/NA49 Collaboration |^,|0| and E802/E859/E866 |l06|,|TT|,|lT| analyze 



in the nucleon-nucleon center of mass frame which is similar (but for asymmetric 
collision systems not identical) to the LCMS at mid-rapidity. 

4.1.5 Square roots of 2, 3, and 5 

Confusion can easily arise when comparing fitted HBT radius parameters with 
the rms or hard sphere radii of the colliding nuclei. As is well known, the 3- 
dimensional rms radius -Rrms,3d and the hard sphere radius -Rbox = 1-2 A^'^ fm 
are related by a factor ^3/5: 

-^rms,3d = {^ ) = {^o + ^s + ^l) = pR. Z = 7-^box- (67) 

The HBT radii are 1- dimensional rms radii (e.g. B^ = (x^)) and thus another 
factor \/3 smaller than -Rrms,3d- If a cold spherical nucleus in its ground state 
could be induced to emit pion pairs, one would thus measure 

Tjcold ^rms,3d ^box ,r,Q\ 

We will call the ratio of the actually measured Rs to this naive expectation the 
"expansion factor" ^. 

At high energies it is often useful to compare with the 2-dimensional rms radius 
of the nuclear overlap region in the transverse plane, -Rrms,2d = (^^o+^s)^ ) which 
is a factor v2 larger than the corresponding sideward radius Rs = -R_l. 

^.2 A Measured Correlation Function 

Collision systems with light projectiles (S-l-S and S-|-Pb at the SPS and Si-|-Au 
at the AGS) have been studied by several experiments, with some systematic 
differences in the results. NA35 included a factor of 2 in the exponent of the 
Gaussian fit function, which yields R parameters larger by v2 than those from 
other experiments. 

Figure n shows three-dimensional correlation functions of Tr"*" and K'^ pairs in 



14.6 A GeV Si-|-Au collisions, measured by E859 [111|. Shown are the projections 
in the three Cartesian directions, with narrow cuts on the momentum differences 
in the other two directions. However, the fits are performed in three dimensions, 
not upon the projections. The solid line shows a Gaussian fit in the CM system 
(y=1.25) according to ( ^3| ) without the cross-terms; the shape of the data is 
reproduced quite well. 

The fit parameters extracted by E859 from pion and kaon correlations are 
shown in Table |^. The radius parameters from kaons are considerably smaller 
than those from pions, but the A parameter is larger |106| , p-lCl . Both trends 
are expected from the fact that kaons are less affected by resonance decays than 



pions [112|; from the discussion in Sec. 2.4.1, however, the effect on the radius 



parameters should have been weaker. It was also postulated that, due to their 



smaller cross section with nucleons, kaons may freeze out earlier than pions [113| 
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Figure 6: vr+Tr^ and K^K^ correlation functions for 14.6 A GeV Si-|-Au colli- 
sions from E859 |111]. The analysis was done in the participant center of mass 
(y=1.25). For the projections the remaining two components of the relative mo- 
mentum were cut to 5 MeV/c < (?perp < 35 (65) MeV/c for pions (kaons). The 
solid lines show a 3-d Gaussian fit. (Figure taken with permission from [|111||.) 



reflecting a smaller source if the latter expands. The E859 analysis [111| showed, 
however, within a fit of reduced dimensionality (2 instead of 3 radius parame- 
ters) that also the pion correlation radii decrease systematically with increasing 
transverse mass of the pair. However, the interpretation of such a 2-dimensional 
fit is not straightforward. 



Jf-.S Asymmetric Collision System,s: First Signs of M±- Dependence 

The first 3-dimensional and liT-dependent correlation analysis was achieved in 
sulphur-induced collisions at the SPS. Table |2| summarizes the R parameters 
extracted by NA44 from a 3-dimensional analysis of pion and kaon correlations 
in 450 GeV p+Pb and 200 A GeV S+Pb collisions at the SPS [|o|,|ll. They 
were extracted from a fit to 

C{q,K) = 1 + A(K)exp(-i?2(K)g2 _ Rl{K)ql - Rf{K)qf) . (69) 

One sees the same trends as observed by E802/859 at lower beam energy |110| , 
111]: the R parameters for kaons are consistently smaller than for pions and the 
correlation strength A is larger. That the R parameters are larger in S-|-Pb than 
in p+Pb collisions should be expected. However, Rg in S-|-Pb collisions is also 
(much) larger than the projectile, even for kaons: for ^^S we should compare 
Rs to Rhox/V^ = 3.8fm/2.23 = 1.7 fm. Thus the system must have expanded 
significantly before freezeout. 

This is supported by the additional observation (see Figure [^ that all three R 
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pair 


Ri 


Rs 


Ro 


A 


K+K+ 


2.75±0.15 
1.71±0.14 


2.95±0.19 
2.09±0.20 


3.77±0.13 
2.07±0.16 


0.65±0.02 
0.83±0.08 



Table 1: E859 correlation function fit parameters [|111[| . Pions have rapidities 
1.2 < y < 1.8 and transverse momenta 100 MeV/c < p± < 800 MeV/c, while 
kaons are accepted for 1.0 < y < 1.7 and 100 MeV/c < p± < 900 MeV/c. 







Rs 


Ro 


Ri 


A 


S + Pb 


7r+ 


4.15±0.27 


4.02±0.14 


4.73±0.26 


0.56±0.02 


P + Pb 


7r+ 


2.00±0.25 


1.92±0.13 


2.34±0.36 


0.41±0.02 


S + Pb 


K+ 


2.55±0.20 


2.77±0.12 


3.02±0.20 


0.82±0.04 


P + Pb 


K+ 


1.22±0.76 


1.53±0.17 


2.40±0.30 


0.70±0.07 



Table 2: HBT parameters in the LCMS measured by NA44 for 450 GeV p+Pb 



and 200 A GeV S+Pb colhsions [[L04| , |114[ . Pions are measured in 3.2 < y < 4.2, 
kaons in 2.7 < y < 3.3. The p± acceptance is < p_L < 0.6 GeV/c for both. 



parameters show a strong dependence on the transverse mass M± of the pairs. 
In 1 88 1 they were compared to a common l/\/M±-law (see Figure |^) which si- 
multaneously reproduces the two pion points and one kaon point in each of the 
three Cartesian directions. This points to a dependence on M± (rather than on 
mo and K± separately), as expected for a thermalized expanding source p9i . 

NA44 compared correlations of K~^ and K~ pairs, to investigate possible effects 
of their different cross sections with nucleons on the freeze-out distribution. No 
significant differences were found, indicating that at SPS energies the KN cross 
section plays a subdominant role in the freeze-out process. This is to be expected 
from the small measured nucleon/pion ratio: p/tt~^ = 0.12 [115|; Ktt scattering 
dominates and is similar for K^ and K~ . 

The fit function (|69[) did not include the cross-term Roi- It is expected to be 



small near mid-rapidity where the NA44 acceptance is concentrated. Subsequent 



analyses showed [89| that including the cross-term changes the R parameters by 
< 10%, less than the 20% systematic uncertainty on their absolute values. 

NA35 has measured the rapidity and transverse momentum dependence of the 
three source radius parameters Rg, Ro, and Ri for central S-l-S, S-|-Cu, S-|-Ag, 
and S-I-Au collisions [105|, by fitting the correlation functions to Eq. ( p9| ) in the 
nucleon-nucleon center of mass. The rapidity dependence of the transverse radius 
parameters is minimal. Ri depends strongly on the pair rapidity and is well 
described by Ri ~ l/cosh(y— y^^,^), indicating approximately boost-invariant 
longitudinal expansion ||6^ as was predicted for collisions at these energies [116 



If this y-dependence in the CM frame were exact, Ri would be independent of Y 
in the LCMS. One should remember, however, that the analysis of |105|| was done 
without including the Roi cross-term which was shown |47| to become sizeable in 
the CM at large Y -Y^^. 

All systems show a similar dependence of Ri on the transverse pair momentum 
K±, with Ri decreasing with increasing K± at a similar rate as found by NA44 
(see above). This provides further evidence for longitudinal expansion of the 
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Figure 7: M_|_ -dependence of the three Cartesian radius parameters in 200 A 



GeV S+Pb colhsions measured by NA44 [B8|. The open (filled) symbols indicate 
values extracted from tt^tt^ {K^K~^) correlations. The solid line is given by 
R = 2 im/y/MZ (M_L in GeV). 



source (see Sec. 2^), and this behaviour is also predicted by hydro dynamical 
simulations of the collisions [92, 105]. 

Compared to Ri, the transverse parameters Kg, Ro from the NA35 analysis 
|105| ] show a weaker i^_|_-dependence, indicating that the longitudal and transverse 
expansions differ. Such a tendency is generically expected from hydrodynamic 



simulations of the collision dynami cs [P2 ,93| and from the model studies presented 
in Sec. 2^. As discussed in section |3.1.1| , the apparently stronger ir_|_-dependence 
of the transverse R parameters observed by NA44 in S-|-Pb collisions |^^ can be 
understood in terms of the different experimental acceptances. 



4-4 Au+Au Collisions at the AGS 

Collective transverse expansion of the emitting source was also observed at the 
AGS. A characteristic K±^ dependence of the transverse radius parameter R±^ was 
found in Au+Au collisions at 11 GeV/nucleon by E866 [gg]. Decreasing values of 
all three Cartesian source parameters with increasing K±^ are observed in Au + 
Au collisions at energies as low as 2 GeV/nucleon by the E895 collaboration [p9|| , 
although at their lowest beam energy of 2 Gev/nucleon the K±^ dependence of 
Rg and Ro seems to disappear. Though the hadron densities at such energies are 
much less than at 158 GeV/nucleon, they are still quite large. Even here, the 
hadronic scattering can generate pressure and cause the source to expand. In 
fact, pair momentum dependence of the fit parameters was already observed in 



streamer chamber data at the Bevalac by Beavis et al. [117| in a 1-dimensional 
analysis of pion correlations from 1.8 GeV/nucleon Ar+Pb collisions. Already at 
that time the dependence was interpreted as evidence for expansion of the source 



before freeze-out [117|. 
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4.5 Pb+Pb Collisions at the SPS 

4.5.1 Cartesian parametrization 

Figure ^ shows projections of the 3-diniensional correlation functions for vr^ and 
vr"*" pairs from 158 GeV/nucleon Pb+Pb cohisions, measured by NA44 (closed 
points) [p9|. The left column shows a 1-dimensional analysis in ginv for com- 




50 100 150 50 100 150 50 100 150 50 100 150 

Qi„„(MeV/c) QTs(MeV/c) 0^0 (MeV/c) Ql (MeV/c) 



Figure 8: Comparison of NA44 Pb+Pb vr^vr^ correlation data (solid circles) 
and RQMD predictions (open triangles). The 1-dimensional projections of the 
3-dimensional correlation function are averaged over the lowest 20 MeV/c in the 
other momentum differences. (Figure taken from 



parison with older data from other colliding systems. Two angular settings of 
the spectrometer cover p_L-ranges of 0-0.4 GeV/c and 0.3-0.8 GeV/c. The fitted 
source parameters, with and without the Rqi cross-term, are given in Table |89|. 



As in the case of sulphur-induced collisions, all R parameters are seen to become 
significantly smaller as K±^ increases, again pointing to collective longitudinal and 
transverse fiow of the source. At small i^_|_ Kg is once again much larger than the 
corresponding value (3.1 fm) of a cold Pb nucleus, indicating that the transverse 
flow also leads to transverse growth, of the collision zone before freeze-out. 

The cross-term i?^^ is non-zero for all data sets, and it is rather large for the 
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pair {{K±)) 


A 


Ro (fm) 


Rs (fm) 


Ri (fm) 


Ri, (fm^) 


iT-7r-{^ 170) 


0.495±0.023 


4.88±0.21 


4.45±0.32 


6.03±0.35 




TT+TT+i^ 170) 


0.569±0.035 


5.50±0.26 


5.87±0.58 


6.58±0.48 




7r+7r+(Ri 480) 


0.679±0.034 


4.39±0.18 


4.39±0.31 


3.96±0.23 




TT-7r-{^ 170) 


0.524±0.026 


5.35±0.25 


5.07±0.35 


6.68±0.39 


10.7±2.9 


7r+7r+(« 170) 


0.658±0.035 


5.98±0.23 


6.94±0.48 


7.39±0.40 


28.1±3.5 


7r+7r+(Rs 480) 


0.693±0.037 


4.59±0.21 


4.71±0.36 


4.15±0.25 


3.1±1.4 



Table 3: Cartesian radius parameters from a Gaussian fit to NA44 correlation 
functions for Pb+Pb collisions [g^ using the Coulomb wave correction. The fitted 
results with and without the cross-term R'^, are shown. {{K±) in MeV/c.) 



low-i^_i_ 7r+ pairs. As explained above, this term should be non-zero in the LCMS 
frame (where these data were analyzed) if the source is not reflection symmetric 
in beam direction. Since the NA44 low-K± acceptance is slightly forward of 
midrapidity, this condition of symmetry is not fulfilled. The fitted R parameters 
all become larger when the cross-term is included. The higher K± acceptance is 
nearer midrapidity, and the cross-term is indeed smaller. 

It should be noted that the fit parameters for positively and negatively charged 
pions differ in an apparently significant way. However, calculation of the x^ P^r 
degree of freedom between the two measured correlation functions yields a value 
of 450/440 |8^. As this is near unity, the experimenters concluded that the 
correlation functions do not, in fact, differ. This illustrates an important system- 
atic limitation in extracting source parameters from Gaussian fits to measured 
correlation functions. Such problems are certainly exacerbated when compar- 
ing data from different experiments where statistical and systematic errors de- 
pend differently on q. These limitations hold regardless of the choice of source 
parametrization, but a comparison of different parametrizations (in which the 
available q-space is differently populated) via the cross-check relations given in 
Sec. 2^ could provide an estimate for the corresponding systematic uncertainties. 

Figure M also shows pion correlation functions calculated from the RQMD event 



generator |86|, using the same charged multiplicity as selected by the experiment 
and a filter simulating the NA44 acceptance. RQMD predicts source size param- 



eters which are slightly larger than the measured ones [B9|, but agrees remarkably 
well with the general trend of the data. It reproduces the larger values of Ri com- 
pared to the transverse R parameters for low K± and predicts radius parameters 
similar to the measured ones at high K±. It overpredicts, however, significantly 
the value of the correlation strength A. This discrepancy is likely due to non- 
Gaussian distortions of the correlation function by pions from resonance decays 



and the resulting systematic uncertainties in Gaussian fits (see Sec. 2.4.1 ). 



4.5.2 Yano-Koonin-Podgoretskii parametrization 

An analysis of charged particle correlations from 158 A GeV Pb+Pb collisions 
with the Yano-Koonin-Podgoretskii parametrization (|30|) (amended by a correla- 
tion strength parameter A) was performed by the NA49 Collaboration [107 , 118|. 
The analysis is not based on identified pions, and the non-pion contamination 
(which contributes to the mixed-pair background but not to the correlated pairs) 



34 



ULRICH HEINZ & BARBARA V. JACAK 



reduces the value of A significantly. However, the values of the YKP fit parame- 
ters are affected only at the 2-6% level [ [107 1. 

The left part of Figure ^ shows the fit results for the parameters i?|| , R± and Rq, 
both as a function of pair rapidity Y = I^tt for small transverse pair momentum 
0.1 < K± < 0.2 GeV/c and as a function of K± for forward moving pairs at 
3.9 < y < 4.4. Rii and R± peak at midrapidity, the maximum of Rn (~ 7 — 8 



0.1<K^< 0.2GeV/c 



3.9<Y^^<4.4 





0.4 0.6 

Kj^ (GeV/c) 

Figure 9: Left: YKP radius parameters for pion pairs from 158 A GeV Pb-|-Pb 
collisions measured by NA49, as functions of the transverse momentum K± and 
rapidity Y = 1^,^ as indicated in the Figure. Right: Effective source rapidity Y^^^ 
for pious as a function of pair rapidity Y = Yj^t^, for pairs with small transverse 
momenta (0.1 < K± < 0.2 GeV). The dashed horizontal line indicates the ex- 
pectation from a non-expanding source. Filled circles are measured data, open 
circles are reflected about midrapidity {Yt^.,^ = 2.9). (Figure taken with permission 
from jlOTl .) 

fm) being somewhat larger than that of R± (~ 6 fm, more than twice as large as 
expected from the initial nuclear overlap region). Since at midrapidity Ri = Rn, 
the NA49 YKP radius parameters RiuR± can be compared with Ri,Rs from 
NA44; they are consistent. While all three YKP radius parameters are decreasing 
functions of K±, the i^^-dependence of R\i is clearly stronger than that of R±, 
indicating dominant longitudinal and somewhat weaker transverse expansion. 

The right diagram in Figure ^ shows a strong correlation between the effec- 
tive source rapidity Y^^ extracted from the YKP fit and the pair rapidity y^j^- 



Following the discussion in Sec. |2.3.1| this is again evidence for very strong longi- 
tudinal expansion of the source. The data points in the Figure seem to indicate a 
slope slightly below unity, as expected from thermal smearing effects at low K± 
(see Figure P) . 

The data were compared [107] to the expanding source model (|4^). In each Y- 



bin the K^-dependence of the i?-parameters can be successfully described by the 
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model, and in particular the strong i^_|_-dependence is reproduced very well by the 
assumption of boost-invariant longitudinal flow. The rapidity dependence of R±, 
however, cannot be reproduced by a constant transverse Gaussian radius R of the 
source; it requires R (and thereby the average transverse flow) to decrease away 
from midrapidity [^,|6^. Near midrapidity one finds for the model parameter 
8 fm, shrinking to i? ss 7 fm at 3.9 < F < 4.4 
3.7±1.6GeV-i 



,10^ ■ From (g^) 
107 1 . For freeze-out temperatures 



Rin {]^ R 

one obtains the ratio rj'j/T 

in the range 100-140 MeV this implies transverse flow rapidities rjf of 0.6-0.72, 

corresponding to average transverse flow velocities of 0.5-0. 6c. 



4.5.3 Emission duration 



In sulphur-induced collisions at the SPS |88,105| and in Pb+Pb collisions ana- 
lyzed by NA44 with the Cartesian parameterization [^ the emission duration 
was found to be very short - consistent with 0-2 fm/c (with considerable sta- 
tistical and systematic uncertainties). NA49 found from their YKP fit to the 
Pb+Pb data a non-zero emission duration of approximately 3 fm/c |107]. All 



these numbers are rather short compared to the emission duration predicted by 
RQMD, namely 3-6 fm/c for S-hPb Q and 7-8.5 fm/c for Pb-hPb @, despite 
the fact that RQMD provides a good representation of the particle distributions. 
This illustrates the difficulty and model dependence of extracting the emission 
duration in relativistic heavy ion collisions, discussed in the previous sections. 

This is different for low-energy heavy-ion collisions at E/A=30-80 MeV where 
emission durations of up to 1400 fm/c were measured |119(1 - the typical evap- 
oration time of a compound nucleus. Such large times can be extracted with 
relatively much less model uncertainty. If the creation of a quark-gluon plasma 
led to a very long-lived intermediate stage near the critical temperature Tc for 
hadronization, it could emit hadrons from the surface over much longer periods of 
time than presently measured. This might leave more easily interpretable traces 



in R^^Q ov R^ [ggll. 



5 COMBINING SINGLE- AND TWO-PARTICLE SPECTRA 



It has been shown that the single-particle m_|_-distributions of pions and heav- 



ier hadrons reflect a transversely expanding source |65, 72, 120|. Such spectra 



are sensitive to a different combination of T and rjf [70, 72|. Combining this 
information with that from the two-particle correlation functions thus allows a 
separation of the collective and thermal momentum components described by ijj 
andT p|,|63|, pll , p^ . 

NA49 has analyzed negative hadron and deuteron spectra in Pb-|-Pb collisions 
at the SPS and combined the fit parameters with their correlation function anal- 
ysis near midrapidity 



|107] ]. Figure 10 shows the allowed regions of freeze-out 
temperature and radial velocity as bands of ztu around each of the three inde- 
pendent fits. These favor a narrowly defined overlap region with T = 120 it 12 
MeV and rjf = 0.55 ib 0.12 [107]. From the model parameter, ijf, one can calcu- 



late the average transverse flow velocity {v±) (see below), given the density and 
velocity profiles. The fit to the h~ spectrum in |l07t| was made without explicit 
consideration of the contributions from resonance decays and heavier hadrons; 
once these are included |12S], the band labelled "/i~" in Figure 10 no longer 
bends over at small T, and the crossing region is shifted slightly downward. 
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Figure 10: Allowed regions of freeze-out temperature T and transverse flow ra- 
pidity 7]f (here denoted as /3_l) for central Pb-|-Pb collisions near niidrapidity. 
The bands labelled h~ and d indicate fits to the negative hadron and deuteron 
spectra; the third band stems from a fit to the sideward radius parameter from 
2-pion correlations. (Figure taken with permission from [[107|.) 



95- 



[ 1241 gave T 

Pb-|-Pb h~ spectra and correlations |125| yielded T 



A simultaneous fit to all the available one- and two-particle distributions was 
done on Si-|-Au data from the AGS by Chapman and Nix |6^]. They found a 
freeze-out temperature T ~ 90—95 MeV and an average radial flow velocity of 
{v±) ~ 0.34c. A comprehensive analysis of single-particle spectra from Au-|-Au 
collisions at the AGS (for which the 2-particle correlations still await publication) 
confirms [64| the low freeze-out temperature T ~ 90 MeV, with an even larger 
average transverse flow {v±) ~ 0.45 c. A preliminary simultaneous analysis of 
spectra and correlations from Pb-|-Pb collisions at the SPS measured by NA44 
100 MeV and (7;_l)_^0.34c. A simultaneous fit to the NA49 

100 MeV with r// = 0.6, 
corresponding to an average transverse expansion velocity {v±) ~ 0.5 c. These 
freeze-out temperatures are somewhat lower (and the radial expansion velocity 
correspondingly higher) than the values extracted from single particle spectra 
alone p|, pOtp^ . 

Knowing the average transverse flow velocity {v±) and the transverse size R 
at freeze-out, a lower estimate for the total expansion time can be obtained. 
Comparing for the Pb+Pb collisions at the SPS the 1-d rms radius at midrapidity, 
i? ~ 8 fm (see Sec. 4.5.2), with the corresponding value of a cold Pb nucleus, 3.2 
fm, one arrives with {v±) < 0.5 c at a lower limit of 10 fm/c for the total duration of 
the transverse expansion. This is larger than the value tq ~ 8 fm/c extracted [jlOTl 
from the measured midrapidity value of R\\ using (^). This illustrates the already 
mentioned problems of interpreting tq directly as the total expansion time, and 
it also indicates that the longitudinal expansion was not always boost-invariant, 
but that the system underwent longitudinal acceleration before freeze-out. 



6 GLOBAL TRENDS 
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6.1 From p+p to Pb+Pb Collisions 

Analysis of correlation functions in p+p and n/K+p collisions yields source pa- 



rameters R of approximately 1 fm [127,128|. Proton collisions upon nuclei yield 



larger effective sources. As shown above, Rg from p+Pb collisions at 450 GeV is 
approximately 2 fm, smaller than the corresponding value (3.2 fm) for a cold Pb 
nucleus. The density of produced particles in such collisions is not large and so 
significant expansion should not be expected; the increase relative to pp is more 
likely due to cascading of the struck target nucleons. 

In heavy-ion collisions at both AGS and SPS energies significant flow is ob- 
served via the M± dependence of the source parameters. Consequently, Rs mea- 
sured with pairs of finite M± provides only a lower limit to the actual size of the 
source. Even so, the measured values reflect effective sources which are already 
considerably larger than the initial nuclear overlap region. For Si+Au collisions 
at 14.6 GeV/nucleon, the ratio ^ of Rs to Rf^'^ (as defined after Eq. (|68|)) shows 
that the source expands by at least a factor ^ = 1.8. As their measured correla- 
tion functions agreed with those from the RMQD model, the E814 collaboration 
used RQMD to correct for the flow effects |^^. They inferred an expansion fac- 
tor ^ = 2.7 by comparing Rt = \/{x'^ + y'^) from RQMD with the 2-dimensional 
radius of the projectile in the transverse plane (see Sec. 4.1.5|) . 



For Au+Au collisions at the AGS, the lower limit on the expansion factor 
from Rg is ^> 1.6, while for Pb+Pb at CERN it is ^ ~ 2. Data from symmetric 



collisions of Fe, Nb, and La at the Bevalac |129, 130(| at 1.3-1.7 GeV/nucleon 
indicate expansion factors S, of at least 1.2-1.4. 

6.2 Beam Energy Dependence 

Central Au+Au collisions at 2, 4, 6, and 8 GeV/A were measured by E895 [^] 



and analyzed with the parametrizations described in Sec. 2.2.2. Though the 
analysis is still preliminary, some very interesting trends are evident. The K±- 
dependence of Rs characteristic of radial flow is observed in all but the lowest 
energy collisions, where pion statistics limit the measurement. Longitudinal flow 
(via R\\{K±)) is seen at all energies. The transverse expansion factor ^ ranges 
from 1.5 to 2.0. It is intriguing that the largest expansion is at the lowest energy, 
where the relative importance of the radial flow appears to be smallest. It is 
tempting to conclude that lower expansion velocity coupled with large final size 
implies a long-lived source. However, the data indicate that the emission duration 
in all cases is quite short, leaving one to wonder why surface emission of pions 
appears to be missing. 

Direct comparisons with Bevalac data are complicated by the inability to select 
the most central collisions. However, the fits use the function C{q±,qH,q^) = 
l + Xeyip[{-q'j_R'j_-qJR?,-{q^)'^T'^)/2], so -R_l may be compared with Rs. Studies 



of 1.7 GeV/A Fe+Fe, 1.5 GeV/A Nb+Nb |T2g] and 1.3 GeV/A La+La collisions 
|130| ] yield R± values (corrected by \/2 for comparison with the AGS fits) of 2.8, 
3.4, and 3.2 fm, respectively. These give ^ = 1.2-1.4. The disagreement with the 
low-energy AGS results may be due to the centrality difference of the collisions. 
Collecting all of the data together, one may look for trends with ^/s in the 
region between 2 and 20 GeV covered by present data. For asymmetric heavy ion 
collisions with small projectiles, Rs or R± is always 3-4 fm, regardless of ^/s. For 
small symmetric systems {A < 100), Rs is likewise approximately independent 
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of ^/s, in the range 2.5-3.5 fm. For symmetric systems with A > 100, it is 
more instructive to look at the expansion factor ^ which appears to increase with 



y/s. This conclusion rehes heavily, however, upon the La+La measurement [130| 
which is low (^ = 1.2). 

7 CONCLUSIONS AND FUTURE PERSPECTIVES 

7.1 Where Do We Stand? 

In this review we have described the development of a sophisticated framework 
with which to extract physics from two-particle correlation measurements in 
heavy-ion collisions. Theoretical and experimental progress in the past decade 
allows characterization of the particle source in these rather complex systems and 
gives access, for the first time, to the dynamical evolution. At last the promise of 
elucidating the space-time evolution of the particle source directly from measured 
quantities has been realized. These new techniques are now also being applied to 
e~^e~ collisions. 

Two-particle correlations measure collective flow of the matter, via the pair 
momentum dependence of the homogeneity region, thus fixing the ratio of the 
freezeout temperature and average flow velocities. Combining this information 
with an analysis of the single particle spectra uniquely separates temperature 
and flow. The detailed characterization of the final state which is now possible 
provides stringent constraints on models simulating the dynamical evolution of 
the reaction zone. 

7.1.1 Heavy-ion collision dynamics 

Analysis of correlation functions has shown that tremendous expansion of the 
system takes place before the hadrons decouple. Longitudinally the source at 
freeze-out features approximately boost-invariant flow while the transverse dy- 
namics is slightly weaker. Still, the transverse radius of the particle emitting 
source approximately doubles from initial impact to freeze-out. These features 
indicate action of a significant pressure, though the hadronic observables are not 
able to indicate which degrees of freedom are responsible for its build-up. 

Already at Bevalac energies, below 2 GeV/nucleon, the hadronic matter ex- 
pands. However, both the radial flow velocity of the matter and the final source 
size tend to increase with increasing beam energy. High energy collisions at the 
SPS result in hadron sources which develop average transverse flow velocities 
of 0.5 c; the hadrons freeze out at temperatures near 100 MeV. Kaons and pions 
flow together, and the observations are consistent with freeze-out from a common 
source. From the onset of expansion to freeze-out, the transverse radius increases 
by a factor of 2.5; given the difference in expansion velocities, the longitudinal 
growth should be about twice that. The total source volume in Pb-|-Pb collisions 
thus grows by a factor of 30! 

7.1.2 Initial conditions 

The single and two-particle distributions are consistent with formation of a ther- 
malized, flowing hadron gas. The behavior of such a gas may be used to ex- 
trapolate back to the early times in the collision, using the measured freeze-out 
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conditions and flow gradients constrained by the data. The freeze-out temper- 
ature, via the equation of state of an ideal resonance gas, provides an estimate 
for the local energy density at decoupling (~ 80 MeV/fm^ at T ~ 100 MeV), the 
measured transverse flow velocity provides the Lorentz contraction factor 7^ to 
correct for the kinetic energy of the expanding matter in the lab frame (7^ w 1.3 
for {v±) ~ 0.5). This gives a freeze-out energy density of about 100 MeV/fm^. 
With the total expansion factor 30 given above, the initial energy density, at the 
onset of transverse expansion, should have been of the order of 3 GeV/fm'^, i.e. 
well above the critical energy density e^ < 1 GeV/fm^ for color deconfinement as 
given by lattice calculations [|131| . Since the initialization of transverse expansion 
requires pressure, this large energy density must have been at least partially ther- 
malized. This estimate of the initial energy density (which is averaged over the 
transverse area of the source) agrees in order of magnitude with simple estimates 



which use the Bjorken formula |116] with the measured multiplicity density and 
assume an initial thermalization time Teq = 1 fm/c. Compared to the latter it has, 
however, the advantage that it replaces several features of the highly idealized 
Bjorken expansion model |116f| (e.g. the parameter Tgq) by measured quantities, 
extracted from the HBT analysis. 

7.2 The Future 

Experimental measurements in the coming years will develop in several directions. 
Let us mention a few in which serious activities are already now visible: 

7.2.1 Heavy-ion collisions at higher energies 

Higher energy collisions (factor of 10 increase in ^/s) will be available at the 
Relativistic Heavy Ion Collider (RHIC) in late 1999. Experiments will use the 
techniques described here to map the freeze-out conditions. Quantifying the 
transverse and longitudinal flow velocities will allow determination of the pres- 
sure build-up in the early stages of the collisions. Armed with the extracted 
velocities, one will be able to back-extrapolate from freeze-out to the time of 
hadron formation. Full analysis of kaon and proton correlations both extends the 
accessible K± range and will verify whether these hadrons are emitted from a 
common source with the pions. 

7.2.2 Azimuthally sensitive HBT analysis 

We have only discussed the analysis of azimuthally symmetric sources. Even if im- 
pact parameter b = never occurs, the "central" event ensemble from which the 
correlation function is constructed is azimuthally symmetric since one averages 
over the orientation of the collision plane. On the other hand, reconstructing the 
latter event-by-event opens up an even richer field of activities and phenomena. 
Azimuthal anisotropies of global event features (in particular directed and elliptic 
flow, for reviews see |132, |133(] ) have been studied since the days of the Bevalac, 



and azimuthally sensitive analyses of single-particle spectra from non-central col- 
lisions have recently attracted a lot of attention |134]. The measured elliptic flow 



of pions and protons at the SPS [135| may play an important role in extracting 



the pressure in the early collision stage |136|. The next question is how these 



azimuthal deformations in momentum-space are correlated with corresponding 
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deformations in coordinate-space. The tools for extracting this information from 
azimuthally sensitive HBT analyses were developed in |137-139|. 



Without azimuthal symmetry of the source, the correlator in Gaussian approx- 
imation is characterized by 6 functions of 3 variables, see (|2^). Following the pio- 
neering studies of Voloshin and Cleland |137], it was shown by Wiedemann [138 



that for sources whose transverse geometric and dynamical deformations have a 
dominant quadrupole component this can be reduced to 6 functions of only 2 
variables {K±,Y). The dependence on the angle ^ between K±^ and the impact 
parameter h is made explicit in the Wiedemann parametrization: 

C(q, K) = 1 ± exp[-i?2g2 _ ^2^2 _ ^2^2 _ 2Rlq,qi] 

X exp[— ai cos$(3gQ -|- q^) + 2a\ 'sva.^qoqs\ 

X exp[-Q2 cos(2$)(5o - g^) -h 2a2 sin(2$)gogs] . (70) 

The parameters Rs,Ro, Ri,Roi are the same as in (P7|); they are given by the $- 
averages of the corresponding functions of K=(Y, K±, <I>) in (|25|) and describe the 
homogeneity regions of the azimuthally averaged deformed source. The two new 
parameters ai^2{K±,Y) are related to the first and second harmonic coefficients 
of Rs in the azimuthal angle $. They were shown [|138|| to characterize dynamic 
and geometric elliptic anisotropies in the source, respectively. 

7.2.3 The average phase-space density at freeze-out 

Determination of the phase-space density of particles at freeze-out will be impor- 
tant to check whether the collisions at higher energies produce a source consistent 
with local thermal equilibrium. It has been speculated that the higher densities 
of produced particles might cause a pion condensate j^. Here, a comparison of 
pion and kaon phase-space densities may be particularly illuminating. 

The average phase-space density of pions at freeze-out has already been studied 
at AGS and SPS energies and found to be consistent with particle emission from 



a source in local thermal equilibrium. Using Eq. (59), E877 analyzed positive 
and negative pions from 10.8 A GeV Au-|-Au collisions near beam rapidity. They 
found |2^ that transverse flow is not needed, and the data are well-described 
by a thermal source. Study of a compilation [Q of data from S-|-nucleus and 
Pb-|-Pb collisions at the SPS showed a nearly universal behavior of the phase- 
space density of pions at freeze-out. The data are again consistent with freeze-out 
from a locally equilibrated source. However, the p^-dependence of the spatially 
averaged phase-space density at mid-rapidity is less steep than that of a Bose- 
Einstein distribution. Such a deviation is expected due to transverse flow. 

7.2.4 Three-pion correlations 

Higher order Bose-Einstein correlations will be studied to seek evidence of cohe- 
rence in the source. An analysis as described in Sec. B^ was done by NA44 [140|, 



albeit due to limited statistics only in one dimension. A three-pion correlation 
function C^iqs) = 1 + Xsexpi-R^Ql), with Q3={pi-p2)^+iP2-P3)^+{P3-Pi)'^, 
was constructed from central S-|-Pb collisions at 200 A GeV; it is shown on the 



left side of Figure 11. The upper dashed line shows the three-pion correlation 
function expected if the source is totally chaotic and symmetric, while the lower 
shows the case of vanishing r-^ (see Sec. pB). The right side of the Figure shows 
that the normalized three-pion correlation function r^ (which is experimentally 
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Figure 11: Left panel: Three-pion correlation function from S+Pb collisions mea- 
sured by NA44, as a function of Q3 (see text). Right panel: Half the normalized 
true 3-pion correlator r^ (here called uj) as a function of Q3. (Figure taken with 
permission from |14C[|.) 



determined as a weight factor called cu ]]140[]) deviates from 1 in the experiment 
in the accessible Qs-range; if this behaviour persisted down to Q3=0, this would 
imply partial coherence of the source. It is hoped that the much higher statistics 
expected at RHIC will settle this important question. 
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